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APRESENTACAO 


Estas notas tratam do cálculo integral de funções de varias variáveis. 


Aprendemos no curso de cálculo de uma variável que o conceito de primitiva 
de uma função contínua e positiva está ligado à noção de área da região 
limitada pelo seu gráfico e o eixo Ox. Por definição, a área é igual à integral de 
Riemann da função e a conexão mencionada acima se dá por meio do Teorema 
Fundamental do Cálculo. 


O objetivo deste curso é generalizar a integral de Riemann para várias variáveis 
(duas ou três). No caso de várias variáveis veremos que o conceito de integral 
dupla está relacionado ao volume de uma região do espaço tridimensional 
limitada pelo gráfico de funções. 


A integral tripla é usada para a obtenção do volume e massa de sólidos mais 
gerais. 


Veremos que para efetuar o cálculo da integral dupla ou tripla usa-se um 
procedimento semelhante ao “Princípio de Cavalieri” estudado na Geometria 
Espacial. Para calcular um volume de um determinado sólido, a ideia é subdividi- 
lo em fatias e calcular a área de cada fatia. Em seguida, calcula-se a integral (ou 
“soma”) da área das fatias. Este é, descrevendo de modo bastante simplificado, 
o significado do Teorema de Fubini, que diz que o cálculo da integral dupla (ou 
tripla) se faz por meio da integral repetida. 


Há também outras maneiras de generalizar o cálculo integral. 


A primeira delas é calcular a integral de funções definidas em objetos mais 
gerais no espaço. 


Por exemplo, suponhamos que um pedaço de arame, feito de um material cuja 
densidade não é constante é descrito como uma curva no espaço tridimensional. 
Qual é a massa do objeto? Para resolver este tipo de problema é desenvolvido 
o conceito de integral curvilinea. 


Analogamente, se tivermos uma placa que tem a forma de uma superfície 
não plana, (por exemplo, um pedaço de um cilindro) feita de um material de 
densidade variável, qual é a massa desta placa? 


Aqui a generalização se faz com o desenvolvimento do conceito de integral de 
superfície. Em particular, se considerarmos a densidade igual a um estamos 
calculando a área da superfície. Por exemplo, iremos aprender como encontrar 
a área de uma calota esférica. 


Ao final iremos estudar as várias generalizações do Teorema Fundamental do 
Cálculo em diversos contextos. 
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O primeiro contexto relaciona o calculo da integral curvilinea sobre uma curva 
fechada com o calculo da integral de superficie na regiao delimitada pela 
curva. (Teorema de Green e Teorema de Stokes). 


O segundo relaciona o calculo da integral de superficie em uma superficie 
fechada e limitada com o cálculo da integral tripla na região delimitada pela 
superfície. Teorema de Gauss ou da Divergência. 


Estas notas foram escritas para ser utilizadas em um curso a distância. Para 
isto elas se dividem em aulas. Cada aula se abre com uma lista de objetivos que 
o aluno deve ter em mente ao estudar o material. Alguns exercícios simples, 
presentes no texto, têm o objetivo de testar se o material exposto está sendo 
assimilado. No final de cada aula, uma lista de exercícios mais elaborada e 
complexa tem o objetivo de ajudar o aluno a absorver melhor os conceitos 
expostos e aprender a operar com eles. 


Belo Horizonte, março de 2012. 
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Integral dupla 


AULA 1: INTEGRAL DUPLA 


Objetivos 1.1 Os objetivos desta Aula sao: 
e introduzir o conceito de integral iterada; 
e calcular integrais iteradas em regiões planas limitadas por gráficos; 
e introduzir o conceito de integral dupla; 


e usar o Teorema de Fubini para calcular a integral dupla por meio da integral iterada; 


1.1 Introdução 


Dois problemas fundamentais são estudados no curso de Cálculo em uma variável: 
1. calcular a área de uma figura plana 


2. encontrar uma primitiva para uma função contínua. Isto é: dada uma função 
contínua f : (a,b) > R, encontrar uma função F tal que F'(x) = f(x), Vx € (a,b). 


Um dos mais belos e úteis resultados que se aprende naquele curso é o Teorema Funda- 
mental do Cálculo (ver [4]), que relaciona esses dois problemas. A área da região plana 
limitada pelo gráfico de uma função contínua e positiva f : [a,b] > R, o eixo -Ox e as duas 
retas verticais x = a e x = b é calculada do seguinte modo: Primeiramente encontra-se 
uma primitiva F(x) de f(x), em seguida, calcula-se a diferença F(b) — F(a). Em símbolos: 


b 
f T a NOE 


Com esta ferramenta em mãos obtemos várias aplicações importantes. O objetivo deste 
curso é generalizar estas ideias para várias variáveis. 

Neste capítulo apresentamos os conceitos de Integral Iterada (ou repetida), de integral 
dupla e as técnicas de integração em duas variáveis. 

De modo análogo ao caso de uma variável, veremos que as integrais duplas são úteis em 
várias situações, entre elas citamos: 


1. cálculo da área de figuras planas mais gerais; 
2. cálculo da área de uma superfície; 


3. obtenção da massa de uma placa plana feita de um material que possui densidade 
variável; 


4. cálculo do volume de uma região do espaço tri-dimensional limitada pelo gráfico de 
uma função positiva de duas variáveis, pelo plano xOy e por planos verticais. 


Veremos que para efetuar o cálculo integral em várias variáveis é fundamental o domínio 
das técnicas de integração em uma variável. Recomenda-se, portanto, que seja feita uma 
breve revisão da integral em uma variável e também de geometria analítica. 


AULA 1: INTEGRAL DUPLA 
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1.2 Integral Iterada 


O procedimento utilizado para o cálculo da integral dupla é a integral iterada que, por ser 
conceitualmente mais simples, é o nosso ponto de partida. Antes, porém, vamos recordar 
o conceito de soma de Riemann de uma função real: 

Se f : [a,b] > R é uma função real, então uma partição, denotada por P, do intervalo [a, b] 
é uma escolha de um número finito de pontos no intervalo: a = zo < 121 <... < Zn = b. Em 
cada sub-intervalo [x;, 1,41] escolhemos um ponto arbitrário, vi € [x;,1;+1] e calculamos 
a sua imagem f(x*). A soma de Riemann de f com respeito à partição P é definida por 
s(f,P) = Dão Fama — 24). 

Se denotamos por |P| = maxf|x;+1 — x;|+ o tamanho da partição, então, para uma função 
contínua f temos 


b 
[ Fede = im 56?) 


Este resultado será usado com frequência nas próximas seções. Um ponto a ser destacado 
é que não importa o ponto que escolhemos x? no sub-intervalo [x;, £;+1] para formarmos 
a soma de Riemann, o limite existe e é igual à integral de Riemann. 

Passemos agora ao conceito de integral iterada: Considere uma função contínua f : [a,b] x 
[c,d] > R. Para cada x € [a,b] fixo, definimos a função contínua de uma variável (y), 
f* : [c,d] > R por f*(y) = f(x,y). f? é a restrição de f ao segmento vertical {x} x |c, d]. 
Segue, o que f” é uma função contínua. Assim sendo, podemos calcular a integral 
definida: A(x =f? f*(y)dy para obtermos uma função contínua que depende apenas da 


variável x. PR existe a integral e A(x)dz. 


Definição 1.2 Chama-se Integral Iterada a seguinte integral T { i f* (y)dy}dx 


Exemplo 1.3 Calcule a seguinte integral iterada: 


[iff aran 


Solução: Observe que, como está indicado, primeiramente integramos em relação a y. A 
variável x é, portanto, mantida constante. Logo 


3 9, oy? Sto Sa 19r 
zy dy = E = 
2 


Como foi dito acima, esta é uma função de x. Calculemos agora a integral de + em 


relação a x para obter: 
[ a => / ae = Be _19 
0 6 


Observação 1.4 Interpretação Geométrica da integral iterada: 

Suponhamos que a função f seja positiva no retângulo R = [a, b] x |c, d]. Então, de acordo 
com o que foi visto no Cálculo em Uma Variável, fixado x € [a,b], a função A(x) é igual 
à área da região plana limitada pelo gráfico de f” e o eixo horizontal entre as retas y = c 
e y = d, conforme a figura a seguir. 
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A integral iterada é a integral da função área A(x) entre os extremos z =a e x = b. 
O próximo exercício será útil para que se compreenda bem a idéia de Integral Iterada: 


Exercício 1.5 Calcule a integral iterada fo [ de 3x + 4y?dy]dx 


Procedendo de modo análogo, definimos a integral iterada 


/ É f EE 


Primeiramente fixamos y € [c,d] para obter uma função de x. A integral B(y) = 
f? f” (x)dx define uma função que depende apenas da variável y. Em seguida, integramos 


B(y) para obter um número: fe B(y)dy 


Exemplo 1.6 Calcule a seguinte integral iterada: 


[uf xy” da] dy 


211 


do: r 5 ‘li 2 2 è 
Solução: Primeiramente calculamos f, cy*dx = =| = & Em seguida, calculamos 
0 
pipas = 8 _ 19 
22> A. e. GS 


AULA 1: INTEGRAL DUPLA 
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1.2.1 Integral iterada em regiões mais gerais 


Até agora calculamos a integral iterada em retângulos [a,b] x [c,d]. Porém, procedendo 
da mesma forma como acima, podemos calcular a integral iterada em regiões mais gerais 
do plano. Por exemplo, se g(x) e h(x) são funções contínuas definidas num intervalo [a, b] 
tais que g(x) < h(a) em (a,b), podendo haver igualdade nas extremidades do intervalo, 
então faz sentido calcular a seguinte integral 


h(a) 
f LF (y) dy, 
g(x) 


para obter A(x) uma função apenas da variável x. Em seguida, podemos calcular f° A(x)dx 
para obter um número. Este é o significado da integral iterada 


1 1 / F=(y)dyjaz 


Observe a ordem em que calculamos a integral. Vejamos um exemplo: 


Exemplo 1.7 Calcular a seguinte integral iterada: 


2 


fy ea ial de 


Solução: Observe que no intervalo (0,1) temos x? < x? com igualdade nos pontos extre- 
mos do intervalo. Primeiramente integramos em relação a y: 


x 2 
T 
f vydy = 50’ 
23 z 


Em seguida, calculamos 


Exercício 1.8 Calcule a seguinte integral iterada i, oe + 2y3)dydz. 


Exercício 1.9 Faça um esboço da região plana R descrita pelas desigualdades: 


lI<a<2, 1l-2<y<l+uw 


1.3 A Integral dupla em retângulos 


Considere uma função contínua f : [a,b] x [c,d] > R definida em um retângulo R = 
a,b] x [c,d]. A maneira de definir a integral dupla que adotaremos é inteiramente análoga 
ao caso uni-dimensional. 

Iniciamos com partições a = ao < a1 <..<an = b e c = ĉo < c1 < ... < Cm = d de cada 
um dos intervalos [a,b] e [c,d], para obter uma partição do retângulo R em sub-retángulos 
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Rij = [ai 0444] x [bj, bj+1], para i = 0,...,n — 1 e j = 0,...,m — 1. Esta partição será 
denotada por P. Escolhendo um ponto arbitrário z;; = (x;,y;) € Rij em cada um dos 
sub-retângulos podemos formar a soma dupla: 


n-1m-1 


YY Has) — 05) (bjt b) 


i=0 j=0 


Se denotarmos por Ax; = (ai+1 — a) e Ay; = (bj+1 — b;), então a soma dupla acima pode 
ser escrita na forma: 


n=1m-=1 
P)= Y f(2)Az¡Ay,. 
i=0 j=0 
Seja [P| = max{Az;, Ay;), onde maz significa que estamos tomando o valor máximo 


entre todos os i =0,...,n—1lej=0,...,m-— 1. 
A Integral Dupla de f no retángulo R é definida, como no caso da integral de Riemann na 
reta: 


f [Foda = timp sos FP) 


Observação 1.10 Para entender o significado do símbolo dA, que chamamos o elemento 
de área em coordenadas cartesianas, basta verificar que para a função constante f(x,y) = 1 
a integral dupla é igual à área do retângulo R. Isto porque, para cada partição P, a soma 
s(f,P) é constante e igual à soma das áreas dos sub-retângulos R;;. Sendo assim, temos: 


| fia=(-ata-o) 


Vejamos mais um exemplo : 


Exemplo 1.11 Calcular f Sia bIx[e,d] g(x)h(y)dA para g e h funções contínuas. Para uma 


partição P qualquer dos intervalos [a,b] e [c,d] , se tomarmos pontos (x;,;) no interior 
de cada sub-retângulo, teremos f(x;,y;) = g(x;)h(y;). Logo a soma s(f,P) se escreve: 
n=1m-1 
s(f,P) =D" glz )hlyj) (wens — z:)(uj+1 — ys). 
i=0 j=0 


Esta soma pode se escrita como o produto de duas somas de Riemann em uma variável: 


n-1 m-1 
s(f, P) = bD g(xi)( Ti+1 — Ti) ID AC h(y \(Yj+1 = y;)). 
1=0 j=0 


Tomar partições de tamanho cada vez menores (isto é com |P| — 0) significa introduzir 
mais sub-intervalos cujos comprimentos tendem a zero ou seja, cada fator (1;+1 — £i) e 
(yj+1 — yj) tende a 0. 

Portanto, usando a definição de integral de Riemann em uma variável, encontramos: 


lima E )h(y)dA = timyppsos(f, P) = canta 


AULA 1: INTEGRAL DUPLA 
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Por outro lado, 


if "ai i ale / 1 À (fs f 1 / "g(a)h(y)daldy 


Conclui-se assim que, para o caso particular em que a função f é um produto de duas 
funções contínuas, uma g(x) que depende apenas da variável x e outra h(y) que depende 
apenas da variável y, o cálculo da integral dupla coincide com a integral repetida. Mais 
ainda, não importa a ordem de integração. 


Exemplo 1.12 Vejamos um caso particular do exemplo acima: 


a ia ay 3 1 3 
cuda= | f zydydz= = xy s 
fo 1 41 2 E 2lo 2 2 4 


Será que o fato observado acima pode ser generalizado? Ou seja: é sempre verdade que o 
cálculo da integral dupla de uma função contínua é realizado por meio de uma integral 
repetida, em qualquer ordem? A resposta a esta pergunta é sim! (ao menos para funções 
contínuas). No caso em que a região R é um retângulo é o que nos diz o 


Teorema 1.3.1 ( Fubini) Se f : [a,b] x [c,d] > R é uma função contínua, então 


HH... id ii fy f(x, y)dylda = [of endas 


Temos assim um poderoso instrumento para a calcular integrais duplas: O cálculo da 
integral dupla é feito calculando-se a integral repetida. Mais ainda, podemos escolher a 
ordem de itegração de acordo com a nossa conveniência. 


Exemplo 1.13 Encontre a integral f fpr SdA para R limitada pelas retas: x = 1, 1 = 2, 
y=ley=2. 
Solução: Pelo Teorema de Fubini 


fg (io 
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1.4 A integral dupla em regiões mais gerais 


Suponha agora que R seja uma região compreendida entre o gráfico de duas funções 
contínuas isto é: 


R=ta<r<bp(r)<y< o(a). 


Uma região R deste tipo é chamada região de tipo I. 
Neste caso, para cada partição a = ay < ay < ... < an = b do intervalo [a,b], obtemos uma 
subdivisão do conjunto R em faixas, denotadas por R; limitadas por trechos dos gráficos 


de 61 (1) e de da(x). 


REGIÃO DO TIPO | 


Cada uma dessas faixas pode ser subdividida em sub-regiões R;; limitadas ou por seg- 
mentos horizontais que correspondem a valores y = bj ou pelos gráficos de qi(x) e de 
do(x). De qualquer modo, sabemos calcular a área, denotada por A;;, de cada uma 
dessas regiões R;j. Escolhendo, em seguida, pontos (x;,y;) € Ri; formamos a soma 
s(f,Q) = So Vito f (Ti Yi) Ai Ao denotarmos por |Q| a maior das áreas Ajj, a 
integral dupla de uma função contínua f : R — R é definida por: 


/ f f(x,y) = limyg¡os(f, Q) 


De modo análogo procede-se para regiões entre dois gráficos R = {c < y < d} e qi(y) < 
x < ¢o(y). Uma região R deste tipo é chamada região do tipo II. Primeiramente, 
obtemos uma partição por meio de uma subdivisão da região em faixas horizontais. Em 
seguida, subdividimos as faixas horizontais. 


REGIÃO DO TIPO II 
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Continuando a analogia com o caso unidimensional, podemos generalizar o conceito de 
integral dupla para funções contínuas por partes. 


Definição 1.4.1 Dizemos que f : R > R é continua por partes se 

a) podemos escrever R = Rı U R5U...Ry, como uma união finita de regiões, 

b) as regiões R; são sub-conjuntos do plano limitados por gráficos de funções contínuas; 
c) quando duas regiões intersectam, a intersecção coincide com as suas fronteiras, isto é 
pedaços de curvas; 

d) a restrição de f a cada uma das regiões é contínua. Ou seja, são definida por desigual- 
dades do tipo: 

Wily) < x < paly) e dir) < y < Yala). 


Definição 1.4.2 Definimos a integral dupla de uma função contínua por partes 


f semis f sendas f Honda dot f EA 


Pode-se demonstrar, mas náo o faremos aqui, que para calcular esta integral dupla basta 
usar o Teorema de Fubini em cada uma das sub-regiões R;. Vejamos em um exemplo como 
fazer isto. 


Exemplo 1.14 Calcule a integral 


f [even inte 


onde R é o triângulo de vértices A = (—1, —1), B = (0,0) eC = (1, —1). 

Primeira Solução: Inicialmente observe que os lados do triângulo são dados pelas retas 
y=2,y=—xey=-—l. Para cada y fixado —1 < y < 0 temos que y < x < —y, ou seja, 
temos uma região do tipo II. Para calcular a integral procedemos como na fórmula acima 
utilizando o Teorema de Fubini. 


| [c+ a= [eur nanty= f Eura 


0 3 3 2 
y? — y 2y* ¡0 
L 2 y) dy 2 la 


Segunda Solução: Podemos ver a mesma região R como uma região do tipo I. Para 


isto consideramos a subdivisão do triángulo dado em dois subtriângulos R; e Ro tais que 
R= R, U Ro, onde 


Ry = {(2,y) ER? -1<2<0, -l<y<a} 


ut=yY 
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Ry = {(a,y) ER? O< 2 <1, -1<y<-s). 


Temos então utilizando primeiramente a Definição 1.4.2 que: 


f [e+ da= f | @v+D dA+ | f (v+ da 


A 
a a 
B=(0,0) 
Xx 
A = (-1,-1) C =(1,-1) 
y= =] 


Em seguida, utilizando o Teorema de Fubini temos: 


f [e+= f [evry ato = [ eh E, dx = 
0 g3 x a e 0 5 
f ++) de = (+ +0), E 


= 
f | curvas f [ut dy t= [rias 


Finalmente adicionando os dois resultados anteriores obtemos - + 3 = 1 o que mostra que, 
neste caso, podemos calcular a integral como uma região do tipo I ou como uma região do 
tipo II. Veremos na Seção 1.5 que isto nem sempre acontece. 

1.4.1 Propriedades da Integral Dupla 


O Teorema de Fubini nos dá o caminho a ser seguido para efetuar o cálculo da integral 
dupla: 


e Primeiro passo: descreva a região de integração por meio de inequações. Observe se 
há necessidade de decompor a região em sub-regides mais simples. 


e Segundo passo: escreva a integral repetida e observe se a ordem pode ser importante 
para facilitar o cálculo das primitivas. 


e Terceiro passo: calcule a integral repetida. 


Algumas propriedades são úteis para o cálculo da integral dupla e são similares ao caso 
de uma variável. 


1. Lia + g(x,y)|dA = Fl Fx, ,yidA + J falo, yda, o mesmo valendo para 
a diferenca. 
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2. f fpcflo,ydA=cf fp f(x,y)dA, c € R constante. 
3. Se f(x,y) > O então J Jr fo, y)dA >0 


4. Se R= BUC com BNC ou é vazio ou é uma união de curvas contínuas. Então 


| | tenaa= | [rodas f | tenaa. 


1.5 Mudança na ordem de Integração 


Na seção anterior, vimos como fazer a integração em regiões do Tipo 1, que são descritas 
por desigualdades: 
a<v<b dim) <y < da); 


e também em regiões do Tipo 2 que são descritas por desigualdades: 
exy<d ily) <2 < plz). 


Há regiões que podem ser descritas tanto como regiões do Tipo I, como do Tipo II. Con- 
sidere o seguinte exemplo: 


Exemplo 1.15 Considere a região R interior ao círculo de raio 1 dado por a? + y? = 1. 
Resolvendo para y temos: y = +v1 — «2, Esta região pode ser descrita como uma região 


do Tipo I, onde d1(1) = — v1 — a? e ¢o(x) = V1 — r? : 
-I<a<l —yl-22?<y< y1-z2?. 


Mas resolvendo para x temos x = +,/1 — y?. E então podemos descrever a mesma região 


como uma região to Tipo II para 41 (y) = —\/1—y? e paly) = y1— y? : 
-1 <y<1 eyi=psasy Iii 


Explorando esta dupla maneira de expressar uma região do plano temos uma boa ferra- 
menta para resolver integrais duplas. Observe o seguinte exemplo: 


Exemplo 1.16 Mude a ordem de integração para resolver a seguinte integral: 


1 v1-a2 
/ V(1—y?) dy dz 
o Jo 


A integral iterada acima é equivalente a integral dupla 


[VE dy do 


onde S é o conjunto dos pontos (x,y) do plano zy onde 0<x<1e0<x< y1-— y? ou 
seja S é o o quarto do disco de raio 1 em que x e y são positivos. 
A região S pode também ser descrita assim 


O<y<1 0<x<yl1-yj?. 
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Temos portanto: 


[fC ver dy a= ff" ya de dy= 


1 1-2 1 31 2 
La Vi y= | a-w- 
0 


0 


y 


Vejamos um outro exemplo deste método 


Exemplo 1.17 Mude a ordem de integracáo para resolver a seguinte integral: 


1 pl 
/ f exp(x”) dx dy 
0 Jy 


Observe a figura que descreve a região R onde a integral iterada acima está definida. A 
integral acima é equivalente a integral 


J | ew dy dz, 


onde R está descrita pelas desigualdades 


O<y<l y<zx<l. 


Mudando a ordem de integração temos 
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1.6 EXERCICIOS 


1. Calcule as seguintes integrais iteradas: 


(a i senle) (1 + y)dydz 


) 

) 

(c) ls, es eV dydx 

(d) fo PA + ydydz 
) 
) 


f ie 32 2+Y dyda 


2. Em cada um dos itens abaixo, faca um esboco da regiáo plana R descrita pelas 
desigualdades: 


(a) l<x<2,1=2<y<1l+.2 
(Oo) laa Sleep 
US al 

3. Em cada um dos exercícios abaixo é possível calcular a integral como uma regiáo 
do tipo I ou como uma regiáo do tipo II. Calcule-as das duas maneiras e verifique, 
em cada caso, que o resultado encontrado é o mesmo. Esboce sempre a regiáo 


triangular correspondente e calcule as equações das retas, lados dos triángulos, pois 
caso contrário, é praticamente impossível resolver o problema. 


(a) 
Je + y)dxdy 


onde D é a regiáo do plano xy dada pelo interior do triángulo de vértices 
A = (0,0), B = (—1,1) e C = (—1, —1). 


(b) 
[Le + y)dxdy 


onde D é a região do plano xy dada pelo interior do triângulo de vértices 
A = (0,0), B = (1,1) e C = (1, —1). 

(c) f fpy(= + y)dxrdy onde D é o triângulo do plano xy de vértices A = (—2,0), 
B = 02) C = (20): 


4. Em cada um dos itens abaixo, primeiramente tente calcular a integral repetida. Em 
seguida, faça um esboço da região de integração e escreva a integral iterada na ordem 
inversa. Observe que as integrais obtidas são fáceis de calcular. 


(a) T ip e” dady 
(b) fo (2? + sen(y?)]dyde 
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5. Para testar a sua intuição, calcule f Je f(z,y)dA para o caso em que f(x, y) = g(x) 
só depende de x e o caso em que f(x, y) = h(y) só depende de y. 


Verifique que: 
a) f Spg(a)dA = [2 g(x) de] (d—c). 
b) J Sn h@)dA = [JE h(y)dy](b — a). 


6. Procedendo como no Exemplo 1.11, utilizando a definição de integral dupla, verifique 
que 


b d 
[rre f taa + hoas- a) 
[a,b] x [c,d] a € 
Aplique esse resultado para verificar que: 


y! [sen(x) + y]dA = E E 
(0, ] [0,1] 2 6 


7. Em cada um dos itens abaixo, calcule a integral dupla f f r f(x, y)dA para a função 
e região especificadas: 


(Gy fy) = sente Rea US y <T 

(b) f(x,y) =x? e R=: {0 < x < cos(y);0 < y < $}. 

Y = =: R a região limitada pelas retas y= vv = 1, x = 2 e y = Q. 

(D lr S er para R a região descrita pelas seguintes desigualdades: 
0O<y<le3y<zr<3 


y 
(e) f(x,y) = x? + sen(y?) para R a região descrita pelas seguintes desigualdades: 
v sysile0srs i 
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Mudanca de variaveis e 
aplicações da Integral Dupla 


AULA 2: MUDANÇA DE VARIÁVEIS E APLICAÇÕES DA 
INTEGRAL DUPLA 


Objetivos 2.1 Os objetivos desta Aula são: 
e apresentar a ideia de mudança de variáveis no plano para calcular integrais duplas; 


e usar as coordenadas polares para calcular a integral dupla em regiões que apresentem 
simetria em relação à origem; 


e apresentar algumas aplicações da integral dupla; 

e encontrar a área de regiões planas; 

e encontrar a massa de placas planares; 

e encontrar o volume de sólidos limitados por gráficos; 


e encontrar a área de regiões contidas em superfícies no espaço; 


2.1 Mudança de variáveis em integral dupla 


O objetivo desta seção é generalizar a fórmula de mudança de variáveis para integrais 
duplas. No caso de uma variável, sabemos que esta fórmula decorre da regra da cadeia. 
Por exemplo, para calcular f sen (x?) x dx, fazemos a substituição de variáveis u = x? de 
modo que du = 2x dx e 


[oo (12) z de = 5 [sen (u) du= > cos(u) = Z cos(a’). 


Quando fazemos uma substituição de variáveis do tipo x = g(u) a integral E f(x)dx 

transforma-se na integral [1º f(g(u))g' (u)du, para g(ua) = ae g(up) = b. Isto é: k f(x)dx = 
Ub / 

[e FO) (de. 

Vejamos, o significado desta fórmula, sem detalhar o argumento: Suponha que g(u) seja 

uma função crescente. Então a uma partição a = zo < 11 < ... < Tn = b corresponde 

uma partição Ua = ug < Uy < ... < Un = Up, para zi = glui), i = 0,1,2....,n. Assim, 

Zi+1 — Zi = g(ui+1) — g(u;). O Teorema do Valor Médio nos diz que existe um ponto 

u; € (ui, um) tal que g(u1) — glui) = g'(už)(ui+ı — ui). Dessa forma, se z; = g(u;), 

então a soma E f(zi)(£i+1 — Zi) se escreve: 


n—1 n-1 
Y Helen- zi] = Y san g(wis1) — gus) = Y F) (ut) (wigs — us) 
1=0 1=0 


Observe que ao tomarmos partições cada vez menores de modo que se Ax; = (x;41— Li) > 
0 então Au; = (uj41 — ui) — 0. Vemos então que o lado esquerdo tende a J? fade 
enquanto o lado direito tende a S F f(g(u))g'(u)du. Em particular, se f(x) = 1, Vx, então 


f dz = Saeg '(u)du, o que justifica a notação dx = g'(u)du 
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Cabe então a pergunta: será que existe uma fórmula semelhante para a integral dupla? 
Ou seja ao fazermos uma mudança de variáveis x = g(u,v) e y = h(u, v), como relacionar 
os elementos de área, dxdy e dudv? 


2.1.1 O determinante como área 


Sejam v = (v1, v2) e w = (w1, w2) dois vetores do plano. No curso de Geometria Analítica 
(ver [5], [6]) aprendemos que a área do paralelogramo de vértices 

O = (0,0), A = (v1, v2), B = (wi, w2) e C = (vı +1, va +w3) é o módulo do determinante 
da matriz cujas colunas são respectivamente, (v1, v2) e (w1, w2): 


v WI 
E 
Ou seja área [OABC] = |ujw2 — wiry]. 
Este resultado pode ser interpretado no contexto de mudança de variáveis, novamente 
usando o que aprendemos no curso de Geometria Analítica. Considere a transformação 
linear (ou mudança de variáveis) do plano definida por (x,y) = T(u,v) = (au+ 8v, yu+dv) 
cuja representação na forma matricial é obtida calculando-se T(1,0) = (a,y) e T(0,1) = 


(6, 9): 
al 


Isto significa que o quadrado de vértices (0,0), (1,0), (0,1) e (1,1) de área igual a 1 tem 
imagem um paralelogramo cujos vértices são, respectivamente, O = (0,0),4 = (a, y), 
B = (8,0) e C = (a + 8, y + ô). 


De acordo com o que acabamos de recordar, a área do paralelogramo OABC 6 igual ao 
módulo do determinante da matriz que representa a transformação ou seja 


área[OABC] = |að — yB| 


Vemos assim que imagem de um pequeno retângulo de lados Au e Av pela transformação 
linear T (u,v) é um paralelogramo gerado pelos vetores Au(a, y) e Av(8, ô). 

Portanto, ao fazermos uma mudança de coordenadas linear, o elemento de área é multipli- 
cado por um fator igual ao módulo do determinante da matriz que representa a mudança 
de coordenadas. Ou seja, se T(R) = R e (x,y) = T (u,v) = (au + Bv, yu + dv), então 


f [ tenaviy= f f Faut poyu+ svlas = dudo 
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Os próximos exemplos ilustram a utilidade da mudança de coordenadas: 


Exemplo 2.2 Calcular f fg rydxdy na região plana R limitada pelas retas x + y = —1, 
+y=1,r—-y=0ex-—y=1. 


Observe que se fizermos u = z + y e v = x — y então nas coordenadas (u,v) a região de 
integração R fica descrita de maneira muitos simples: R’ =: {—1 < u < 1;0 <v < 1}. Ou 
seja R’ é um retângulo, ao passo que R é um paralelogramo. 

Para calcular f f pvydrdy, qualquer que seja a ordem escolhida para calcular a integral 
repetida, deveremos subdividir R em três partes e realizar três integrais. Faça um esboço 
da região para se convencer disto. Observe que se fizermos u = x + y e v = x — y então 
nas coordenadas (u,v) a região de integração R fica descrita de maneira muito simples: 
R=: {—1 < u < 1;0 <v < 1}. Ou seja R' é um retângulo, ao passo que R é um 
paralelogramo. 

A mudança de coordenadas u = xz + y e v = x — y é descrita por uma transformação 
T(u, v) = (442, 52), cuja representação matricial é dada por: 


| 


cujo determinante tem módulo igual a | — ; — 


Portanto, 
Lu v u vil 
f |ru f Ne o 5) o dudu 


[No] [al So] Le 
| Nile 
La 


Exercício 2.3 


Calcular f Je zydxdy na região plana R limitada pela retas 2x + 3y = 0, 2z + 3y = 1, 
—21+y=0e-2x+y=1. 


2.1.2 Integral dupla em Coordenadas Polares 


As coordenadas polares são muito úteis para descrever regiões planas que possuem alguma 
simetria em relação a um ponto (origem). 

Portanto, para usar integrais duplas nessas coordenadas, é importante que se obtenha o 
elemento de área em coordenadas polares. 

Recordando, escolhidos uma origem no plano cartesiano (0, 0) e um semi-eixo ( o semi-eixo 
Ox positivo), um ponto P £ O do plano pode ser localizado por meio de dois números: 
r = dist(P,O) a distância de P a O, e 6 = ângulo orientado positivamente entre o vetor 
OP e o semi-eixo Ox. 
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Segue da definição que se P = (x,y) em coordenadas cartesianas, então r = 12 + y? 
ie > | 

e cos(0) = NESTE Obtemos assim uma mudança de coordenadas definida por (1,0) > 

(x,y) dada por x = rcos(@) e y = rsen (0). 


Pergunta-se então, qual é o elemento de área em coordenadas polares? 


Observe que um pequeno retângulo de área drd corresponde a um pequeno setor circular 
em que dr corresponde a uma pequena variação do raio e d0 corresponde a uma pequena 
variação do ângulo, conforme a figura a seguir 


0 


Responder a esta pergunta equivale a obter uma aproximação para a área deste setor. 
“ a 2 P : 2 

Ora, sabemos que um setor circular de ângulo d e raio r tem área igual a A(r) = 500. 

Portanto, derivando em relação a r, obtemos que o elemento de área em coordenadas 

polares é igual a dA = rdrdé. 

Uma alternativa, usando diferencial, é a seguinte: x = rcos(0), y = rsen(0), implica que 

dx = cos(0)dr — rsen(0)d0 e dy = sen(0)dr + rcos(0)d60. Ou seja dx e dy são obtidos por 


uma transformação linear cuja representação matricial é 


dx \ _ | cos(0) —rsen(0) dr 
dy ) | sen(0) rcos(@) do 
Essa matriz tem módulo do determinante igual a r. Assim analogamente ao caso linear, 


temos que o elemento de área é igual a rdra(0). 
Vejamos alguns exemplos: 


Exemplo 2.4 Talvez o exemplo mais simples seja utilizar coordenadas polares para ob- 
ter a conhecida fórmula da área de um disco D de raio a com centro na origem. Em 
coordenadas polares o disco se escreve r = a. Temos 


f | aa- Fr rárdo = 29| gah 


No próximo exemplo, veremos como calcular a área de uma região plana limitada por 
curvas escritas em coordenadas polares: 


Se C é uma curva plana fechada, contendo a origem no seu interior e definida em coorde- 
nadas polares por r = f(0), para a < 0 < 8 então, a área da região plana limitada por C 


é igual a 
B pico) By á 
/ / rdrd0 -| af (9) do. 
a 0 a 
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Exemplo 2.5 Seja C dada por C =: r = 3 + 2cos(0), 0 < @ < 27. 
A área da região limitada por C é igual a 


27 p(3+2cos(0)) 27 2 
/ / rdrd0 = / (3 + 2cos(9))" ¿y 
0 0 0 


2 
Ou seja 
27 p(3+2cos(0)) 27 1 
/ Il irdi = | 134 2cos(9)2do 
o Jo o 2 
27 1 z 1 27 
f 518 + 2c0s(0)Pd0 = 5 f [9 + 12cos(0) + 4cos?(0)]d8 
0 0 
Recordando que f cos?(0)d0 = f 5[1 + cos(20)]d0 = 5[0 + 5sen(20)], 
obtemos 


T 


27 1 2 
> f [9 + 12cos(0) + 4cos*(0)]d0 = 3198 + 12sen(@) + 20 + sen(20)] ; 
0 
Substituindo os valores 0 = 0 e 0 = 27, 


i a 1 
3 / [9 + 12cos(0) + 4cos*(0)]d0 = 3187 +47] = 11r 


2.1.3 Fórmula da mudança de variáveis em Integral dupla 


Nas seções anteriores, vimos dois exemplos de como se transforma o elemento de área 
quando fazemos uma mudança de coordenadas; o caso linear e o caso das coordenadas 
polares. Neste segundo caso usamos uma aproximação linear das coordenadas. 

Nosso objetivo é mostrar uma fórmula geral para qualquer mudança de varáveis: Se 
(x,y) = F(u,v) é uma mudança de coordenadas da forma x = g(u,v) e y = h(u,v) e 
R' é a região do plano de coordenadas (u,v) tal que F(R’) = R, se gu, gv, hu € hy 
denotam as derivadas parciais, então 


f [ senavay = | f Fotu) hugdu- go oa dudo 


Exemplo 2.6 A fórmula quer dizer que se fizermos a substituição de variáveis x = 


glu, v) = u? — v?, y = h(u,v) = uv então para escrever a integral J [r f(=, y)dxdy nas 


novas coordenadas (u,v), primeiramente calculamos gu = 2u, gy = —20, hy =v € hy =u 
de modo que |guhv — guhy| = 2u? + 2v2. Assim, a fórmula de mudança de variáveis se 
escreve: 


f [ evda = II. f(g(u,v), h(u, v))2(u2 + v?)dudv 


Por que é valida a formula? 

A obtenção rigorosa desta fórmula foge dos objetivos destas notas. Podemos apenas dar 
uma idéia do tipo de aproximação utilizado. 

A ideia é usar a diferencial das funções g e h para fazer uma aproximação linear: dg = 
gudu + gydv e dh = hydu + hydv. 
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Um pequeno retângulo de lados du e dv é enviado pela diferencial em um paralelogramo 
gerado pelos vetores du(gu, hu) e du(gy, hy). 
Mais especificamente, as imagens dos vértices A = (u,v),B = (u + du, v), C = (u,v + dv) 
e D = (u + du,v + dv) de um pequeno retângulo são 44 = F(u,v), By = F(u + du, v), 
Cı = F(u,v + dv) e Dy = F(u + du, uv + du). 


Usando a diferencial: 

Bı = B' = F(u,v) + du(gu(u, v), hulu, v)), 

Cı ~ C = F(u,v) + du(gu(u, v), hy(u, v)) 

Dı x D' = F(u,v) + du(gu(u,v), hulu, v) + du(gy(u, v), hy(u, v)). 


Não entraremos em detalhes sobre o significado desta aproximação (=) que, na verdade, 
está contida na ideia de derivada como aproximação linear da função. 

Intuitivamente, isto significa que estamos substituindo um pequeno trecho da curva que 
passa pelo ponto A’ = F(u, v) pelo seu vetor tangente neste ponto. De modo que a região 
do plano limitada por curvas com vértices 41,B1j,C4 e Dı tem área aproximadamente 
igual à área do paralelogramo de vértices A’, B’,C’ e D’. 


Como vimos, no caso linear, o paralelogramo tem área igual a |guhv — gyhy|dudv. 

Para ter uma melhor ideia do que está acontecendo, examine dois casos particulares para 
a mudança de coordenadas: 

o primeiro F(u,v) = (u,g(u,v)), a mudança de coordenadas leva reta vertical (u = 
constante) em reta vertical e o elemento de área é uma região compreendida entre dois 
gráficos de função u +> g(u,v). Neste caso a afirmação é que a área de uma pequena região 
(isto é, Au pequeno) é aproximada por |g,|dudv. 


Analogamente, se considerarmos mudanças da forma F(u,v) = (f(u,v),v), em que as 
retas horizontais (v = constante) são preservadas. O elemento de área é igual a |fuldudv. 
Veja a figura 
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2.2 Aplicações da Integral dupla 


2.2.1 Área de figuras planas 


Sobre esta aplicação, não há muito o que acrescentar, pois esta foi a motivação para o 
conceito de integral dupla. Para calcular a área de uma região limitada por curvas precisa- 
mos primeiramente descrevê-la como uma região ( ou uma união de regiões) compreendida 
entre dois gráficos e em seguida, calcular a integral dupla por meio da integral repetida. 

Exercícios: Use a integral dupla para encontrar a área das seguintes regiões planas R C R?: 


1. R=: a região limitada pelas retas x = —1, y = —1 e pelas curvas x = (x + i) e 


y=r- r’. 


2. R=: a região limitada pela curva y = 1 — z? e pelas retas y = 2x e y = 0. 


2.2.2 Volume de sólidos limitados por gráficos de funções 


Também aqui a questão é usar a interpretação geométrica da integral dupla: Se f(x,y) 
é uma função positiva, então f f p f(x, y)dA é igual ao volume do sólido M limitado pelo 
gráfico de f e o plano Oy, sobre uma região plana R. 

Interpreta-se assim porque cada parcela da soma 


n m 


s(f,Q)= >> ay) Ay 


i=0 j=0 


é igual ao volume do bloco retangular Bj; de base [a;, a;+1] x [b;,b;+1] e altura f((x;, ys) 
para Aij = (Qi41 = ay) (544 = b;). 

Deste modo, a soma s(f, O) é igual ao volume de um sólido B = UBi;; constituído pela 
união de todos os blocos. 

Veja a figura 


Como, ao tomarmos partições cada vez mais finas, isto é, tais que |Q| é cada vez menor, 
o valor de s(f, Q) converge para f J, f(x, y)dA faz todo o sentido definir 


vol(M) = I f TA 
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Mais geralmente se f(x,y) e g(x,y) são funções contínuas tais que f(x,y) > g(x,y) sobre 
uma região plana R, então, o volume do sólido M limitado pelos gráficos de f(x,y) e 


g(x,y) é igual a vol(M) = f Jalf(a,y) — g(x, y)JdA. 


Exemplo 2.7 Encontre o volume do sólido M limitado pelos planos x = 0, y = 0, z = 0 
ez=1l1-2-y 

Solução: o sólido está compreendido entre os gráficos de z = 0 e z = 1 — x — y e está 
situado no primeiro octante. Portanto, o volume do sólido é igual a f {,(1 — x — y)dA 
para R a projeção ortogonal do sólido sobre o plano xOy. R corresponde exatamente aos 
pontos (x,y) € R? do plano tais que existe um valor z para o qual o ponto (x, y, z) pertence 
a M. A região R é limtada pelas retas x = 0, y = 0 e pela intersecção de z = 1 — £ — y 
com z = 0. Ou seja, a reta 1 — x — y = 0 ou xz + y = 1. Portanto R é o sub-conjunto do 
plano caracterizado pelas seguintes desigualdades: x > 0, y>Oex+y< 1. Eo volume 


do sólido é igual a 
f [0->-ya 
R 


Escolhendo integrar primeiramente em relação a y e depois em relação a x teremos: Fixado 
x, a variação de y, dentro da região R é descrita por 0 < y < 1 — x enquanto 0 < x < 1. 
De modo que 


ffu-a-mas= f f t-e- f wav- pr as 


[ [n-2-vaas f u x —x(1—2) 51 ajda 


Integrando: 


1 1 1 1 1 1 
TJ [1 — z — y]dA = [z — ha? pp +(x) = 
k 2 2 


2 3 o 6 


Exemplo 2.8 Encontre o volume V do sólido situado abaixo do paraboloide z = 4-1? —y? 


a acima do plano z = 0. Temos que: 


v= | f 4-2- ara 
D 


Em que região D do plano devemos fazer a integral? Para entender isto observe que x 
e y variam no domínio D interior a interseção de z = 4 — 2? — y? com o plano z = 0. 
Substituindo obtemos 0 = 4 — x? — y? ou z? + y? = 2, ou seja D é o interior do disco de 
raio 2. 

Para efetuar o cálculo utilizamos coordenadas polares x = r cos 0 y = rsen 0. Obtemos que 
4 = z? — y? = 4 — r°. Logo, 


27 2 
V= TJ 2 — 2? — y drdy = / | (4 — r2)rdrd0 = 
D o Jo 


27 8 32 
= 27(8 — =) = — rT. 
; ml 3) 37 
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2.2.3 Massa de placas planas 


A integral dupla permite-nos calcular massas de placas planas (ou lâminas) não ho- 
mogêneas. O termo placa significa que estamos desprezando a espessura do objeto. 


A densidade é uma propriedade do material e, por definição é a razão entre massa e área. 
Estamos supondo que é conhecida a função densidade p(x, y), positiva e contínua. 


Novamente, para uma subdivisão de uma placa R em pequenas sub-regiões Rij de área 
A;j, escolhido um ponto (x;,y;) € Rij, então o produto p(x;,y;) Ai; é uma aproximação 
para a massa de Rj; . 


Portanto, a soma m(p, Q) = 54 0 p(x, y¡) Aj é uma aproximação para a massa da 
placa. 


Novamente tomando partições cada vez mais finas vemos que a soma m(p, Q) converge 
para f Sr p(x, y)dA que é igual à massa total da placa. 


Exemplo 2.9 Calcular a massa de uma placa plana circular dada por z? + y? < 1, cuja 


densidade é dada por p(x, y) = y 2? + y?. 


Solução: usando coordenadas polares, a placa é dada por 0 < r < 1 e a densidade é 
p(r,0) =r. 


Portanto a massa é igual a Pu i, r2drdo. Calculando a integral, obtemos: 


27 al 27 re 
/ | Para = | [—]od0 
o Jo o 3 
27 1 
f f r2drdo = E 
o Jo 3 


2.2.4 Area de superficies parametrizadas 


Os conceitos desenvolvidos nesta subseção serão utilizados posteriormente quando tratar- 
mos a integral de superfície. No momento, o objetivo é calcular a área de regiões contidas 
em superfícies. Esta é uma generalização importante da área de figuras planas e trata do 
conceito de área de uma região contida em uma superfície parametrizada S. 


O exemplo mais simples de tal superfície ó o gráfico de uma função f(x,y), ou seja o 
conjunto S =: f(x,y, f(x, y))}. Observe que S pode descrita por uma parametrização, que 
a cada ponto (x,y) € R? associa o ponto (x,y, f(x, y)). 


Gráficos de funções de duas variáveis foram estudados no Cálculo III. Vimos que se f é 
diferenciável, então para cada ponto (xo, yo) de seu domínio a equação do plano tangente 
ao gráfico de f no ponto (xo, yo, f((xo,y0)) é dada por 


z= f(20,Yo) + falto, Yo)(@ — xo) + fy(Xo, yo) y — yo). 
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Recordamos como obter esta equação: Por definição, o plano tangente é o plano que 
passa pelo ponto (29, yo, zo) e é gerado pelos vetores Xy, = (1,0, fe(xo,y0)) e Xy = 
(0,1, fy(xo,yo)). No Cálculo Diferencial de várias variáveis, vimos que o produto veto- 
rial Xz x Xy é um vetor normal ao plano. Portanto, a equação vetorial do plano tangente 
é: 


[(z, y, 2) — (Zo, Yo, 20)].[Xz x Xy] = 0 
Ou 
z= f(20,Yo) + fe(Xo, Yo)(x — xo) + fylzo, yo y — yo). 
O nosso objetivo é mostrar como usar a integral dupla para formular o conceito de área 


de uma superfície e calcular algumas áreas. Por exemplo, a área de uma esfera de raio 1. 
Comecamos generalizando a ideia de gráfico como uma superfície. Temos a seguinte 


Definição 2.2.1 Uma superfície parametrizada regular é uma aplicação r : U C R? > RÌ, 
r(u,v) = (x(u, v), y(u, v), (u, v)) tal que os vetores 


rulu, v) = (welt, v), Yu (un, Zulu, v)) e toluço) = (£u (uy luv), 2 (uo) 
sao linearmente independentes para todo (u,v) € U. 


Em outras palavras, o vetor ru(u,v) não é um múltiplo do vetor r,(u,v) ou, equiva- 
lentemente, ru(u,v) x ry(u,v) £ 0. Se na definição acima fixamos uma das variáveis, por 
exemplo fazendo v = vo a imagem de r(u, vo) é uma curva contida na superficie S. Observe 
o exemplo a seguir. 


Exemplo 2.10 Considere a superfície dada por r(u,v) = (2cos(u), 2sen (u), v). Então 
r,(u,v) = (—2sen (u),2 cos(u),0) e ry(u,v) = (0,0,1). Observe que se (x,y,z) são as 
coordenadas da imagem de r(u, v) temos x? + y? = 4 e a coordenada z é qualquer. Temos 
portanto que a superfície dada por r(u, v) é um cilindro. A superfície é regular em todos os 


i j k 
pontos, pois temos ru(u,v) x rv(u,v) =| —2senu 2cosu 0 |= (2cos(u), 2sen (u), 0) Æ 
0 0 1 
0 para todo u. Observe que se fixamos v, por exemplo, v = 4 temos que r(u,4) é um 
círculo contido no cilindro. Se fixamos u, por exemplo u = 7 a imagem de r(4,v) é uma 


reta, que é uma geratriz do cilindro. Veja a figura. 
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Exemplo 2.11 O próximo exemplo é o gráfico de uma função: r(u,v) = (u,v, f(u, v)). 
Temos ru(u,v) = (1,0, fulu,v)) e ry(u,v) = (0,1, fo(u,v)) que são linearmente indepen- 
dentes pois 


1 J 
ru (u, v) x ry(u, v) =i 1 0 falu, v) E (—fu(u,v), —fu(u, v), 1) #0 
01 


é um vetor normal ao plano tangente a superfície no ponto (u,v). 


Exemplo 2.12 Sejar(u,v) = (e“cos(v),e“sen(v),u), então ru(u,v) = (e“cos(v), e“sen(v), 1) 
e ry(u,v) = (—e“sen(v), e“cos(v), 0). 

Se ru(u,v) e ry(u, v) fossem linearmente dependentes, então olhando para a última coor- 
denada, teríamos um absurdo, 1 = 0. 

Logo r(u,v) = (e“cos(v),e“sen(v),u) é uma superfície parametrizada regular. 


Exemplo 2.13 Se r(u,v) = (cos(u)cos(v), cos(u)sen(v), sen(u)), então ru(u,v) = 
(—sen(u)cos(v), —sen(u)sen(v), cos(u)) ry(u,v) = (—cos(u)sen(v), cos(u)cos(v), 0). 


E ralu, v) x ry(u, v) = cos(u)(cos(u)cos(v), —cos(u)sen(v), —sen(u)) 4 0 
para u Æ 5 eux ar. 


Exemplo 2.14 r(u,v) = (ucos(v), usen(v),0), u 4 O é uma superficie parametrizada 
pois, r,(u, v) = (cos(v), sen(v),0) er, (u,v) = (—usen(v), ucos(v),0) e rulu, v) xr, (u, v) = 
(0,0, u) £ 0. 


A condição ru(u,v) x ry(u,v) £ 0 para uma superfície parametrizada regular nos permite 
proceder analogamente ao caso dos gráficos e definir o plano tangente a S em um ponto 
Po = r(uo, vo): o plano que passa em Po e é gerado pelos vetores r, (uy, vo) € Ty (UN, vo). 
A definição de plano tangente nos dá imediatamente que o vetor ry,(uo, vo) X Lo(uo, vo) é 
um vetor normal ao plano tangente. 


Portanto um vetor normal ao plano tangente à superfície no ponto Pp é dado por: 


Ty X Ty = dim = ZuYv), (Elo = ie (Ei = Pa) 


Assim, se designamos xo = x(uo, vo), Yo = Y(uo, vo) € zo = Z(uo, vo), a equação do plano 
tangente se escreve: 


[(x, y, 2) — (xo, yo, 20)]-[Pu(uo, vo) x ru(uo, vo)] = O 


Chea Eni Zu Yw (£ xo) (Luo Luv) y Yo) (Gate Tulio = zo) =0 


Observação 2.15 O plano tangente contém os vetores tangentes a todas as curvas conti- 
das em S e que passam pelo ponto Po. De fato, seja t — (u(t), v(t)) uma curva diferenciável 
tal que u(0) = uo e v(0) = vo. 
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Então, y(t) = r(u(t),v(t)) = (x(u(t), v(t)), y(u(t), v(t)), z(u(t), v(t))) é uma curva no 
espaço tal que para t = 0 temos r(uo, vo) = Po. O cálculo do vetor tangente a esta curva, 
usando a regra da cadeia nos fornece: y'(t) = (u (Dx +0 (t)av, U (t)yu +0 (0)Y,, (0) zu + 
v'(t)z,). Que pode ser escrito na forma Y(t) = u (O) (Lu, Yu, Zu) + Y (E) (Lv, Yv, zv) isto é: 
Y(t) = (u(t)ru(u,v) + v (t)r,(u,v). Fazendo t = 0 obtemos y (0) = u(O)ru(uo, vo) + 
u'(t)ry(uo, vo). Esta expressão significa precisamente que o vetor y'(0) é uma combinação 
linear dos vetores r, (uy, vo) e ro(ug,v)o que são os geradores do plano tangente. Neste 
sentido é que dizemos que o plano tangente é uma aproximação linear da superfície S no 
ponto Po. 


Observação 2.16 Os vetores ru(uo, vo) e ro(uo, vo) geram um paralelogramo no plano 
tangente com vértices Po, Po + ru(uo,vo) Po + Lo(uo,vo) e Po + ru(uo, vo) + Pv(uo, vo). 
Como vimos anteriormente, este paralelogramo possui área igual a |r,,(uo, vo) x Py (Uo, Vo)|- 


Definição 2.2.2 A área de A(D) de uma região D C U contida no domínio de uma su- 
perfície parametrizada regular r(u, v) = (x(u, v), y(u, v), z(u, v)) é igual à seguinte integral 


dupla: 
A(D) = TJ |ru X roļdA 
D 


(dA o elemento de área do plano (u, v)). 


A expressão |ru x rẹ|dA é chamada elemento de área da superfície parametrizada. 
Nos exemplos a seguir, pretendemos mostrar intuitivamente o sentido desta definição e a 
sua motivação: 


Exemplo 2.17 Seja r(u,v) = (au+cv,bu+dv, 0). Calcular a área da região da superficie 
que corresponde à região do domíio (u, v) igual a D = [1,2] x|—1, 1]. Calculamos r,,(u, v) = 
(a,b, 0) e ry (u, v) = (c,d, 0) Portanto |ru x ry| = |ad — be]. 


1 p2 
a(p) = f [ jad = veaa = | f jad — bc|dudv 
D 141 


A(D) = |ad — be|(2 — 1)(1 + 1) = 2/ad — be]. 


Ou seja, 


Exemplo 2.18 Calcular o elemento de área da superficie r(u,v) = (x(u, v), y(u, v), 0). 
Primeiramente calculamos r,(u,v) = (£u, Yu, 9) e ry(u,v) = (tu, Yu, 0). Portanto, |ru x 
ro| = |LuYv — LyYu|. Logo, conforme a definição, 


A(D) = | f leuw- oid 


Observe que esta é precisamente a fórmula da mudança de coordenadas (com f(x,y) = 1) 
já utilizada anteriormente. Isto significa que a noção de área que estamos utilizando é uma 
generalização da fórmula da mudança de coordenadas na integral dupla. Note também 
que estamos usando a mesma idéia intuitiva de que a área do paralelogramo formado pelos 
vetores tangentes a uma superfície parametrizada é o elemento de área da superfície. Esta 
ideia corresponde à ideia de aproximação linear da superfície pelo seu plano tangente. 
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Exemplo 2.19 Encontre a area da esfera de raio 1. 

Solução: primeiramente precisamos de uma parametrização da esfera. Considere a semi- 
circunferência C de raio 1 no plano yOz parametrizada por (0, sen(v), cos(v)), para O < 
v < 7. Ao girarmos a curva C em torno do eixo Oz, obteremos a seguinte superfície 
parametrizada: r(u, v) = (cos(u)sen(v), sen(u)sen(v),cos(v)), 0O<u<2re0<u<r. 


z 


É fácil verificar que os pontos (x(u, v), y(u, v), z(u,v)) que estão na imagem da parame- 
trização satisfazem z(u, v)? +y(u,v)* + 2(u, v)’ = 1 ou seja, a imagem da parametrização 
está contida na esfera de raio 1. 


Também é fácil de verificar que os pontos que não estão na imagem são os pontos do 
meridiano intersecção da esfera com o semi-plano x = 0, y > 0 ( Verifique isso!). 


Assim, usando a simetria da esfera, se calcularmos a área da região D da superfície parame- 
trizada r(u, v) = (cos(u)sen(v), sen(u)sen(v), cos(v)), descrita por0<u< Se0<v<s 
e multiplicarmos por 8 obteremos a área da esfera. 


ry (u,v) = (—sen(u)sen(v), cos(u)sen(v),0) 
ry(u, v) = (cos(u)cos(v), sen(u)cos(v), —sen(v)) 
rulu, v) x ry(u, v) = sen(v)(—cos(u)sen(v), sen(u)sen(v), —cos(v)) 


|rulu, v) x ry(u, v)| = |sen(v)| 4 0 para v 4 0,7. 


A(D) = [ [ sen(v)dudu 


Calculando: 


T 


[ sen(v)dudu = [Prcostonó tas E 


2 
/ [0 + 1]du = = 
0 2 
Portanto, a área da esfera é igual a 85 = 47. 


Exemplo 2.20 Encontre o elemento de área de um gráfico. Seja (x,y, f(x,y)) o gráfico 
de uma função diferenciável f(x,y). 


Vimos que Ely) x ry (5, y) = (—fe(x,y), —fy(z,y), 1) # 0. 
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Logo, calculando a norma deste vetor, obtemos o elemento de area do grafico: 


[rl y) x Ela) = fet + fy? + 1d A. 


Exemplo 2.21 Calcular a área da região do gráfico da função f(x,y) = xy sobre a região 
R limitada pelo círculo da raio 1, x? + y? < 1. 
Como fr = y e fy = x, o elemento de área do gráfico é igual a yy? + x? + 1dA. 


Pelas características da região R e da função que iremos integrar, é conveniente usar 
coordenadas polares. 


Assim, R corresponde à região 0O<r<1,0<0<27. 


Portanto a área é igual a 
27 1 
A(R) = / f Vr? + Irdrdo. 
o Jo 


Fazendo a substituição u = 12 +1, obtemos, du = 2rdr, 0 < r < 1 corresponde al < u < 2 
de modo que 


1 2 
1 1 1 
f V+ irdr= | Vuzdu = zu? |? = =[22 — 1]. 
; or =a 


Logo a área da região é igual a 
27 q 27 
1 2 
/ f Vr2 + Irdrdé = =[22 -1 f do = Z [93 — 1). 
o Jo 3 0 3 


Exemplo 2.22 (Superfície de revolução) Seja C =: (0, f(v), g(v)) uma curva regular 
no plano yOz , isto é, uma curva tal que o vetor tangente (0, f'(v), g'(v)) é não nulo, Vv. 
Se além disso f(u) > 0, então a curva não intersecta o eixo Oz e rotação de C em torno 
do eixo Oz gera a seguinte superfície parametrizada regular (superfície de revolução): 
r(u,v) = (f(v)cos(u), f(v)sen(u), g(v)), O< u < 27. 

O nome "superficie de revolução” (ou de rotação como designam alguns autores), vem do 
fato de que, ao fixarmos v = vo, obtemos a curva r(u, vo) = (f(vo)cos(u), f(vo)sen(u), g(vo)), 
que é um círculo contido no plano z = g(vo), centrado na origem e de raio igual a f(vo). 
Este círculo é chamado paralelo. 

A curva obtida ao fixarmos u = uy é chamada meridiano ou geratriz. Dizemos que uma 
superfície de revolução é gerada pela rotação do meridiano em torno de um eixo (no caso, 
o eixo Oz). 

Verifiquemos o fato da superfície ser regular: 


ru(u, v) = (—f(v)sen(u)), f(v)cos(u), 0) 
ry(u,v) = (f'(v)cos(u), f'(v)sen(u), g'(v)), 


de modo que 


ru (u,v) x rolu, v) = f(v)(9'(v)cos(u), g/(v)sen(u), f’(v)) # 0. 
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Pois |ru X ro] = f(v)1/g/(v)? + f'(v)? é positivo porque f(v) > 0 e g'u} + f(v’? > 0 (a 
curva é regular). 
Concluímos que o elemento de área da superfície de revolução é igual a 


dA = flu yg (w? + f’(v)?dudv. 
Em particular, se g(v) = v, então estamos rodando um gráfico y = f(u) no plano yOz em 
torno do eixo Oz. Neste caso, o elemento de área é igual a 
dA = f(v)\/14+ f(v)dudo. 


Como dA não depende da variável u, podemos integrar uma vez e exibir uma fórmula para 
a área da superfície S obtida pela rotação do gráfico de uma função f(v) em torno de um 
eixo entre dois pontos v =a e v =b. 


=f Piso 14+ f(0) Y + Peau = 20 f 11) TETE 
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2.5 


10. 


11. 


EXERCICIOS 


. Utilize a mudança de coordenadas x = au e y = bu para encontrar a área da região 


plana limitada pela elipse de equação a SF y =l. 


. Encontre o volume do sólido situado abaixo do parabolóide z = 9 — x? — y? e acima 


do plano z = 0. 


. Encontre o volume do sólido situado abaixo do parabolóide z = 2 — x? — y? e acima 


do parabolóide z = x? + y?. 


. Encontre o volume do sólido descrito pelas seguintes desigualdades: x? + y? < 1 e 


r? +y +2 <4. 


. Considere a mudança de variáveis F (u,v) = (x,y) definida por: £ = u, y = u?—v?+v 


(a) Faça um esboço da imagem por F do quadrado [0,1] x [0,1]. 


(b) Calcule o elemento de área. 


. Seja u = z? — y? e v = 2xy uma mudança de coordenadas que envia o quadrado 


R={1<a2<2,1<y< 2} em uma região D' no plano (u,v). 


Calcule f f p dudv diretamente e, em seguida, usando a mudança de coordenadas. 


. Encontre a área da região limitada pela curva dada em coordenadas polares r = 0? 


para 0<0< 5- 


. Se R := {(x,y)/4 < 2? +y? < 9; x > 0;y > 0), use coordenadas polares para calcular 


| iL zydA. 


. Utilize coordenadas polares para calcular 


2 pv2r—x2 
f l (3x — 2y)dydz. 
o Jo 


Use coordenadas polares para calcular f f R eo (+ dA para R(K) o disco centrado 
na origem e de raio K. Observe o que ocorre com o resultado quando K > oo. 


Em cada um dos ítens abaixo, calcular a massa da placa descrita por uma região 
plana R com densidade p. 


(a) R =: {[0,1] x [0,5]} e densidade p(x, y) = 1 + 32? + 5y?. 


(b) R a região limitada pelas curvas y = x e y = 2? e densidade p(x,y) = 47. 
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12. Calcular a área da seguinte superficie: 


r(u,v) = ((2+cos(u))cos(v), (2+cos(u))sen(v), sen(u)) para0d<u<5e0<us<r. 


13. Em cada um dos ítens abaixo, encontre a área da região descrita na superfície: 


(a) R região do plano z = 1 — x — y dentro do cilindro z? + y? = 4. 


(b) R região do parabolóide limitada pelo cilindro x? —22+y? = 0 (use coordenadas 
polares). 


14. Calcule o elemento de área de uma plano dado pela equação ax + by + cz = d. 


15. Encontre a área do tronco de cone obtido pela rotação da reta z = y em torno do 
eixo Oz, paral < y < 2. 


16. Encontre a área da superfície obtida pela rotação da curva (0, em ) entre os pontos 
v = —1 e v = 1 ( catenóide). 
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Integral Tripla 


AULA 3: INTEGRAL TRIPLA 


Objetivos 3.1 Os objetivos desta Aula sao: 


e introduzir o conceito de integral tripla; 


e enunciar o Teorema de Fubini que, analogamente ao caso da Integral Dupla, permite 
calcular a integral tripla por meio da integral repetida; 


e utilizar a integral tripla para encontrar o volume de regiões do espaço limitadas por 
superfíces; 


e usar coordenadas cilíndricas e coordenadas esféricas no espaço R? para calcular a 
integral tripla em regiões com certos tipos de simetria. 


3.1 Integral Tripla em um bloco retangular e o Teorema de 
Fubini 


A definição de integral tripla segue exatamente os passos da integral dupla. As dificuldades 
que geralmente ocorrem estão na descrição ou na decomposição do domínio de integração, 
que é agora uma região D do espaço tri-dimensional cujo bordo é formada por uma união 
de superfícies. 

Repetiremos assim, de modo sumário o que foi feito nas seções anteriores. 

Primeiro passo: a integral tripla em um bloco retangular. Seja D = [a, b] x [c, d] x [p, q] C R? 
um bloco retangular. Seja f : U C R? é uma função contínua definida em um aberto U 
que contém D e P uma partição de D em blocos Dijk = [x;, 2541] X [yj, Yj4+1] X [2k, Zk+1] 
formada a partir de partições de cada um dos intervalos. Ou seja, consideramos tp = 
axa <<. < Tn = b, Yo = c < Yy < Y2: <... <y = dezo =p< z < 
Z2 < ... < Zm = q. Como anteriormente, o tamanho de uma partição P é denotada 
[P| = maz{(zi+1 — i), (Yj+1 — Yj), (Zk+1 — Ze) } € 


s(f,P) =D flat, yj, 2h) (Tiri — ti) (j41 — Yj) (2041 — 2), 
onde (a, y;, 2) é um ponto qualquer do bloco Dix. 


Definição 3.1.1 A integral tripla da função f(x,y,z) no bloco retangular D é igual : 


f | [ tic ar = jim s(t P) 


Note que se f(x,y,z) = 1, então s(f, P) é igual ao volume do bloco D, ou seja (b — a)(d — 


c)(q — p). Logo 
[| fev=e-a@-oa-m, 


o que justifica chamarmos dV de elemento de volume em coordenadas cartesianas. 
Segundo passo: utilizamos o chamado Teorema de Fubini, que de maneira semelhante 
caso bi-dimensional, afirma que para fazer uma integral tripla podemos fazer as integrais 
iteradas que agora serão três. Conforme a ordem em que fazemos a integral repetida 
escrevemos dV = dxdydz ou qualquer uma das permutações de dx, dyedz. 

De maneira mais precisa: 
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Teorema 3.1.1 (Fubini) Se f : U C R? > R é uma função contínua definida em um 
conjunto aberto que contém um bloco retangular D = [a,b] x [c,d] x [p,q] então 


J J [tenaa = ff Hen adote 


De fato, podemos calcular a integral usando qualquer uma das permutações do elemento 
de volume. Entretanto, é preciso atenção para que a ordem em que escrevemos o limites 
de integração seja compatível. 

Assim, no enunciado acima, temos: 


JJ [fewaa = Pf remamiaaz 


Vejamos alguns exemplos de como calcular a integral tripla em um bloco retangular utili- 
zando o Teorema de Fubini: 


Exemplo 3.2 Vamos calcular f f fp f(x,y, z)dV onde D = [0,1]x[0,1]x[0,1], f(x,y, z) = 
xyz. Como a função dada é simétrica por permutações das variáveis, nao há vantagem de 
escolha da ordem de integração. Assim, escrevemos: 


JJ [av= fif if auzarana: 


1 
A primeira integral E gyzda = uate = e. 
1 2,11 
A segunda f “dy =[%F] = 4. 
Final t 1 z dx = z2 i _ 1 
inalmente fo a pS e |= es 


Exemplo 3.3 D = [0,1] x [-1,1] x [1,2], f(z,y,z) = xe”! — zx?. Observe que se es- 
colhermos integrar a primeira parcela na ordem dxdydz, então a primeira integral será 
i xe*Ydx que envolve uma integração por partes. Entretanto, se escolhermos integrar 


7 _e~* bem mais 


A . = 2 1 
primeiramente em relação a y, então teremos: e xe dy = e*Y =e 
-1 
fácil (e rápido). 
Assim sendo escrevemos 


[fo — zz?dV = fui ev! — zx2dyldaldz. 


E claro que já vimos este tipo de escolha na integral dupla. O Teorema de Fubini é im- 
portante porque nos permite trocar a ordem de integração conforme a nossa conveniência. 
Prosseguindo com o exemplo: 


1 
I xe”! — zz’]dy = e” — e” — 22x” 
-1 


1 31 2 
/ e? — e™® — 2zx°]dx = |e” te tS 223) a get SEA 
0 
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E finalmente, 


Cabe agora a pergunta: em que tipo de regiões do espaço além dos blocos retangulares 
podemos definir a integral tripla? 

Vejamos algumas respostas: 

I) O primeiro tipo de região segue imediatamente da integral dupla pois é a região do tipo 
produto R x J, de uma região R plana limitada por dois gráficos de funções e um intervalo 
I = [a,b] 

Analisemos o seguinte exemplo: R =: [(x,y)a < x < b;ġ(x) < y < v(x)), I = [p,q]. O 
cálculo é feito por meio da q repetida (Teorema de Fubini). Observe que fixados x e 
y, então a função F(x a = o É 1 f(x,y, z)dz é uma função contínua. Logo podemos calcular 
a integral dupla f fp F R (x gd como fizemos anteriormente: 


[fr sil [fe reddi 


¡tl 


Exemplo 3.4 Encontre f f f, (1 +y + z)dV, para D definida pelas desigualdades x < 
y<se0<z<1. 

Solugao: Pelo Teorema de Fubini, calculamos a integral iterada em D. Já sabemos que 
0< z< 1. Portanto, basta descrever a região plana R que corresponde aos pontos do 
dominio de integração. Porém, esta região é definida pelas desigualdades x < y < a?. 
Vemos que as curvas y = x ey = x” que limitam a região se intersectam em x = x? ou 
seja para os valores x = 0 e x = 1. Dessa forma, a região D é descrita pelas desigualdades 
O0<z<lae<y<a22?e0<2< 1. Logo, 


fJ [ervas f f Pesos daar: 


Calculamos agora a integral iterada: 


1 pa? pl 1 pa? 21 1 pa? i 
/ f 1 (x+y + z)dzdydz = / f [(x + y)z + =| dydx = / f [x + y + 5]dydz. 
o Je Jo o Jo 2 lo o Jo 2 
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1 px? 1 2 2 
1 z 

[| [ety sldudo = flay +t dx 

ede 2 A a Ola 


1 pa? 1 4 
1 x x 7 
i / Saad i / i" 2 ü lda 30 


II) Para regiões limitadas por gráficos g(x,y) < z < h(x,y) sobre um retângulo R = 


[a,b] x [c,d] o procedimento é semelhante: a integral F(x,y) = no f(x,y, z)dz define 
uma função contínua de modo que o Teorema Fubini implica em 


J | | tenaa f f renia- f fif tensia 


III) Região D do espaço limitada pelo gráfico de duas funções: g(x,y) < z < h(x,y). 
Observe que não está especificada a região do plano que corresponde à integral. Devemos 
primeiramente responder à pergunta: Qual é o conjunto R dos pontos do plano (x,y) que 
correspondem à região D? Melhor dizendo, qual é o conjunto de pontos (x,y) para os 
quais é válida a desigualdade g(x,y) < h(x, y)? 


Exemplo 3.5 Se g(x,y) = z? + y? e h(z,y) = 2— x? — y? e o domínio de integração é 
dado por D := z? +y? < z < 2— a? — y? então estamos interessados em obter o subconjuto 
R := {(2,y)\a? + y? < 2-27 — y?). Ora, determinar R significa descrever os pontos que 
satisfazem à inequação x? + y? + a? + y? < 2 isto é, x? + y? < 1. Ou seja R é o conjunto 
de pontos interiores ao circulo de raio 1. 


No caso geral, procedemos da mesma maneira, primeiramente obtemos o conjunto de 
pontos do plano R que satisfazem à inequação g(x,y) < h(x,y) e, em seguida, usamos 
novamente o Teorema de Fubini para obter: 


J | | siena = f | renea- f fif” tenzia 


Uma recomendação: para fazer a integração neste tipo de região, é útil ter um bom 
esboço dos gráficos das funções envolvidas para a seguir, determinar a região R dos pontos 
do plano que corresponde ao domínio de integração. Observe que R é precisamente a 
imagem do domínio D pela projeção ortogonal m(x, y, z) = (x,y). Logo, é a região plana 
limitada pela curva de equação g(x,y) — h(x, y) = 0. 
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Exemplo 3.6 Encontre f f fpzdV para D a região do espaço limitada pelos gráficos 
x=1- 27%, x = 0, entre os planos y = —1 e y= 1. 

Solução: observe que pela descrição da região de integração D, é mais conveniente escrever 
dV = dxdzdy na integral iterada. 

Assim descrevemos os limites de integração na seguinte forma: fixados y e z a variação 
de x é dada pela desigualdade 0 < x < 1 — z2?. A intersecção dos gráficos é dada por 
1 — 2? = 0 ou seja z = 0 e z = 1. Portanto, as seguintes desigualdades definem os limites 
de integração para a integral iterada: 0 < z <1—27,0<z<le-l<y<l. 


Portanto, 
1 pl q1-22 1 pl 
Jf f-av=f Ll zdzdzdy = | Ja aaa 
D -1J0 Jo -1/0 


1 1 1 22 z4 
f [| ea-2eedy = f E-File = 5 


Exercício 3.7 Como seria a integral repetida escrita na forma dxdydz? 


Roteiro 3.8 Em qualquer uma das situações descritas para calcular a integral tripla, a 
integral repetida é obtida a partir da descrição da região D de integração da seguinte 
forma. Suponha que escolhemos fazer a integração na ordem dydxdz: primeiramente, 
fixamos duas variáveis 7 e z e obtemos uma desigualdade envolvendo gráficos de funções, 
por exemplo, h(z,z) < y < g(x, z). Em seguinda, projeta-se a região D em uma região de 
integração das demais variáveis R. Fixada a variável z, a região R também é descrita por 
uma desigualdade do tipo ¢1(z) < x < ¢2(z). Finalmente, a projeção da região R sobre o 
eixo Oz é um intervalo descrito por uma desigualdade: a < z < b. Temos então: 


b pd2(z) pglz,z) 
| [ [tena ff f Tenant 
D a 1(2) h(x,z) 


3.2 Coordenadas Cilíndricas 


Se usarmos coordenadas polares no plano xOy e coordenada cartesiana usual no eixo 
Oz, obteremos um sistema de coordenadas no R? denominado Coordenadas Cilíndricas 
definido por: 
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O elemento de volume nestas coordenadas é igual produto do elemento área em coordena- 
das polares, isto é dA = rdr dê pelo elemento dz ou seja: 


dV =r dz dr d0 


Quando usamos coordenadas cilíndricas? O critério é exatamente o que usamos para as 
coordenadas polares isto é, quando a descrição da região de integração fica descrita de 
modo mais simples e a integral pode ser calculada facilmente. 

Ao optarmos pelo sistema de coordenadas cilíndricas devemos observar a simetria da 
região de integração em relação ao eixo Oz e a função que iremos integrar ao fazermos a 
substituição de variáveis. 

Os cilindros e os cones são exemplos de algumas superfícies que satisfazem o critério de 
simetria. Vejamos alguns exemplos. 


Exemplo 3.9 O cilindro z? +y? = R?, em coordenadas cilíndricas se escreve r = R. Note 
que a variável z não aparece nesta equação, portanto a figura é realmente um cilindro. 


Exemplo 3.10 O cone z = ay/1? + y? em coordenadas cilíndricas se escreve z = ar. 


Exemplo 3.11 Em coordenadas cilíndricas uma superfície de revolução: 


(f(v) cos(), f (v)sen (0), g(v)), 


se escreve r = f(v), z = g(v), não dependendo da coordenada 6. 


Exemplo 3.12 Utilize coordenadas cilíndricas para calcular o volume do sólido Q limi- 
tado pelos gráficos das funções z = x? + y? e z = 2 — 2? — y?. 

Solução: Substituindo as coordenadas polares nas expressões z = z? +y? e z = 2- x? —y?, 
obtemos z = r? e z = 2—r?. O fato de que as funções não dependem de @ significa simetria 
dos gráficos em relacáo ao eixo Oz. 


2 


A interseção dos gráficos é dada por r? = z = 2—r?. ou seja, 2? = 2, our = 1. 


Observe que o sólido é limitado superiormente por z = 2 — r? e inferiormente por z = r? 


projetando-se no plano xOy no disco r° < 1, conforme a figura abaixo: 


Assim sendo, o sólido Q é descrito pelas desigualdades: r? <z<2-r2,0<r<le 
0<0<2r. 

Portanto, lembrando que o elemento de volume em coordenadas cilíndricas é dV = rdzdrd0, 
temos 
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27 pl p2-r? 
vol(Q) = f / / rdzdrd0 
0 0 Jr2 


27 1 
vol(Q) = / / r(2 — 2r2)drd0 
o Jo 
2T r | 


vol(Q) = / [+= z 10 = 


3.3 Coordenadas Esféricas 


No parágrafo sobre aplicações da integral dupla, calculamos a área da esfera de raio 1 
centrada na origem. Para isso, usamos a seguinte parametrização: 


(9, 4) > (sen (6) cos(0), sen (H)sen (0), cos($)),  para0<0<2r, 0O<g<r. 


Os círculos 0 = ĝo = constante correspondem aos meridianos: (sen(@)cos(09), sen(@)sen(O9), cos(d)). 
As curvas que correspondem a $ = do = constante são círculos de raio igual a sen(do) e 
são chamados paralelos (sen(¢o)cos(@), sen(do) sen(00), cos(do)). 

Estas duas famílias de curvas estabelecem uma maneira de localizar pontos na esfera por 
meio de um par de ángulos: (0, 4), que significa localizar o meridiano e paralelo em que o 
ponto encontra-se. 

Por exemplo, o ponto (3, 3, a) localiza-se no paralelo é = | e no meridiano 0 = Y. 
Analogamente para uma esfera de raio R, de equação cartesiana z? + y? + 22 = R?, ob- 
temos a parametrização (Rsen (¢) cos(0), Rsen (¢)sen (0), Rcos(¢)). Da mesma maneira 
que utilizamos coordenadas polares para obter um sistema de coordenadas no plano, po- 
demos construir um sistema de coordenadas no espaco usando cooredenadas esféricas do 
seguinte modo: fixada uma origem O, se P € R, seja p = dist(P,O). Se P 4 O, então 
p > 0. Nesta situação, o ponto P encontra-se numa esfera centrada na origem e de raio p. 
Para localizar um ponto nesta esfera precisamos encontrar o par de ángulos (0, 4). Dessa 
maneira, podemos atribuir a qualquer ponto P 4 O as coordenadas 


(psen ($) cos(0), psen (@)sen (0), p cos()) 
Exemplo 3.13 Considere o ponto (yv3, 3,2). Obtemos p(P) = 4, 6(P) =% e 0(P) = 3. 


Definimos um sistema de coordenadas esféricas que para cada ponto P € R? associa um 
terno ordenado (p, 0, ¢). 

Como era de se esperar, a equacáo de uma esfera de centro O e raio R nestas coordenadas 
é bastante simples: p = R. 


Exemplo 3.14 Outra superfície de equacáo muito simples em coordenadas esféricas é o 

cone, de equação cartesiana z = ay/ 2? + y?, a > 0. Para obtê-la, basta fazer substituição 

x = psen ($) cos(0), y = p sen (d)sen (0), z = p cos(#) para obter x? + y? = p?sen?(¢) : 
pcos(p) = apsen(¢) 


cot(p) =a 


Equivalentemente q = constante. 
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Para calcular integrais triplas em regides no espaco que envolvam superficies que pos- 
suem algum tipo de simetria em relação à origem, talvez seja vantajoso usar coordenadas 
esféricas. 

Para isso, primeiro apresentamos o elemento de volume em coordenadas esféricas: 


dV = p*sen(@) dp de de 


Esta expressão é obtida de modo similar ao que foi feito com a mudança de coordenadas 
no plano. Não discutiremos aqui a sua obtenção. Apenas lembramos que ela pode ser in- 
terpretada como o volume de um paralelepípedo formado por um paralelogramo tangentea 
esfera e um radial, apontando para fora. Em termos de volume, obtemos o produto do 
elemento de área da esfera pelo elemento dp. 


Exemplo 3.15 Encontre o volume do sólido 2 no interior da esfera de centro na origem 
e de raio 1 e acima do cone z = y £? + y?. 


Solução: Como vimos anteriormente, 2 é limitado por duas superfícies, que em coordena- 
das esféricas possuem equações p = 1 e cot(p) = 1. Ou seja p = 1 e q = F, respectivamente. 
Logo, a região 2 é descrita pelas desigualdades (verifique !): O<p<1,0<0<2re 


0<6<5 
vay f f [ar [7 f f seno ap ae do 
[=P oE ao 
fl JÁ send) dó do = à f costi ad 


va = f f fav- iu- 4 "e i 2 


0 2 


Exemplo 3.16 Encontre o volume do sólido Q limitado pelos cones S1, z = \/x2 + y? e 
So, z = y 3(12 + y?) e pelo plano z = 3. 

Solução: Como foi observado anteriormente, a equação de um cone do tipo z = ay/ 22 + y?, 
a > 0 em coordenadas esféricas é : cot(p) = a. 

Portanto, o sólido Q é limitado por superfícies que em coordenadas esféricas se escrevem: 
Sı := cot(¢) = 1 ou ¢= §, 


So := cot(d) = V3 ou, usando trigonometria, cos(¢) = se isto é, 0 = § 
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e o plano pcos(¢) = 3 (verifique isso!). 


Portanto, as seguintes desigualdades descrevem o domínio de integração: 0 < 6 < 27, 
ps, eves —3_ Logo, o volume é igual à integral repetida: 


cos(p) 


27 T Be 
4 cos(¢) 2 
Vol(Q) = / f f p sentó) dp do do 
o Jz Jo 


Que calculamos 


3 


dia sen (4) dọ do = 


LS 


2T A 9 
E / [ cos3(¢) sen(p) dọ do = 
9 


6 


27 1 =. 


E um bom exercício ( trabalhoso) calcular o volume deste sólido usando coordenadas 


cartesianas. 
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3.4 


10. 


11. 


EXERCICIOS 


. Em cada uma das regiões D abaixo escreva a integral tripla :/ f f p f(£,y,z)dV na 


forma de integral iterada: 


(a) D := |" +974 27 < 1}. 
(b) D a região dentro da esfera x? + y? + z? = 2 e acima do gráfico de z = x? + y?. 


(c) D a regi ao fora do cone z? = 2? + y? e dentro da esfera x? + y? + 2? = 2 


. Calcule o volume do sólido limitado pelos seguintes planos x = 0, y = 0,z=0e 


eo Uap es = Il 


. Use integral tripla para encontre o volume do tetraedro de vértices: A = (0,0,0), 


B = (1,0,0), C = (1,0,1), D = (1,1,0). (Sugestão: encontre a equação do plano 


BCD. 


. Encontre o volume de um tetraedro regular de aresta igual a v2. (Em que pontos 


colocaremos os vértices?) 


. Encontre a seguinte integral f f JoydV para Q = {(2,y,2z)/Va? +y? < z < 1} 
. Encontre o volume do sólido B =[(x2,y,2)/y1?+y2<2<6-=u?*-— y?) 


. Use coordenadas cilíndricas para calcular a integral 


ie a a o Rea zdzdydx 


. Em cada um dos ítens abaixo, esboce a região do espaço R? descrita em coordenadas 


esféricas pela equação: 


. Encontre o volume do sólido limitado inferiormente pela esfera x? + y? + 2? = 4z e 


superiormente pelo cone z = y x? + y?. 


Encontre o volume do sólido contido no octante x > 0, y > 0, z > 0 e limitado pela 
esfera de centro na origem (0, 0,0) e raio 2 e pelos planos verticais x = y e y = V3. 


Encontre a equação do cilindro x? + y? = k? em coordenadas esféricas. 
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AULA 3: INTE 


Integral Curvilinea 


AULA 4: INTEGRAL CURVILINEA 


Objetivos 4.1 Os objetivos desta Aula sao: 
e introduzir o conceito de curva parametrizada; 
e introduzir o conceito de integral curvilínea de uma função escalar sobre uma curva; 
e introduzir o conceito de integral curvilínea de um campo vetorial sobre uma curva; 
e enunciar e demonstrar o Teorema de Green; 


e estudar os campos conservativos; 


4.1 Introdução 


Imagine um pedaço de arame, na forma de uma curva C no espaço tri-dimensional. Su- 
ponha que o arame é feito de um material cuja densidade é uma função que a cada ponto 
p € C, associa um valor f(p). Deseja-se calcular a massa do objeto. A ideia é proce- 
der como foi feito nos capítulos anteriores: usar uma integral para encontrar a massa do 
objeto. Isso será feito, subdividindo C em pequenos pedaços, bem aproximados por seg- 
mentos de retas. Calcula-se a massa de cada pedaço e, em seguida, soma-se para obter 
uma aproximação da massa do objeto. A massa total é o limite dessa aproximações ou 
seja uma integral. 

Nosso objetivo inicial será estender o conceito de integral de Riemann para subconjuntos 
mais gerais: curvas e superfícies no espaço. Neste capítulo faremos a extensão para curvas 
usando os conceitos de caminho e de comprimento de arco e serão tratados os seguintes 
assuntos: 


e caminhos regulares; 

e comprimento de arco de caminhos regulares; 
e a integral em curvas regulares; 

e campo vetorial em R3; 

e integral curvilínea (ou integral de linha); 


e Teorema de Green, que relaciona a integral curvilínea em uma curva fechada plana 
com a integral dupla na região limitada pela curva; 


e campos conservativos no plano. 
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4.2 Curvas Regulares 


Definição 4.2.1 Um caminho contínuo é uma aplicação y : [a,b] > Rê, y(t) = (x(t), y(t), (z(t) 
tal que as funções x(t), y(t) e z(t) são contínuas. Se as funções x(t), y(t) e z(t) são dife- 
renciáveis em (a,b), então dizemos que y é um caminho diferencidvel . 


Nesse caso, o vetor y(t) = (a'(t),y'(t), 2’(t)) é chamado vetor tangente a y no ponto y(t). 
Se y descreve o movimento de uma partícula movendo-se no plano, então o vetor tangente 
é a velocidade da partícula no instante t. 


Exemplo 4.2 ó(t) = (t2,t,0) é um caminho diferenciável com vetor tangente 6’(t) = 
(2t,3t2,0). Por exemplo, para t = 1 temos o vetor 6’(1) = (2,3,0) que é tangente ao 
caminho no ponto ô(1) = (1,1,0). 


Definição 4.2.2 Um caminho y(t) diferenciável é chamado regular, se o vetor tangente 
nunca se anula ou seja, se 7/(t) = (a(t), y(t), 2 (t)) 4 (0,0,0) ou equivalentemente se sua 
( / 


norma não se anula: |y'(t)| = ya Divo ro? £0. 


Definição 4.2.3 Dizemos que C C R? é uma curva regular se é a imagem de um caminho 
regular. 


Em outras palavras C é uma curva regular se existe um caminho regular y : J CR > R? 
tal que y(1) = C. 


Exemplo 4.3 y(t) = (cos(t), sen(t),t) é um caminho regular pois y'(t) = (— sen(t), cos(t), 1) 
e |7/(t)| = ycos?(t) + sen?(t) +1 = V2 


Exemplo 4.4 Por outro lado, o caminho do Exemplo 4.2 , 6(t) = (t?, t?,0) não é regular, 
pois o vetor tangente na origem é nulo. Compare as figuras. 


Y 


x 
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Exercício 4.5 Verifique se o caminho q(t) = (t, t?, t°) é regular ou não. 


Para generalizar a definição de integral de Riemann de uma função f em um caminho, 
primeiramente devemos dar sentido à noção de elemento de comprimento de arco, que 
substituirá o elemento dx na integral simples uni-dimensional. 

Lembre-se que o princípio geral que usamos é: 


1. fazer uma partição P do caminho y em pequenos arcos de comprimento A;s, 


2. em seguida, calcular o valor da função f(t*) em um ponto contido em cada um dos 
sub-arcos e formar a soma: s(f,P) = E FEA 


3. Finalmente tomar o limite da soma quando |P| > 0. 


Vejamos como formalizar um pouco mais essa ideia: 


Definição 4.2.4 Se y(t) é um caminho regular, o comprimento de arco de y entre dois 
pontos é igual a integral 


s(to, t1) = [ y (t)|dt 


0 


Exemplo 4.6 Encontre o comprimento de arco da curva y(t) = (cos(t), sen(t),t) para 
entre os pontos y(0) e 7(4). 

Solução: basta aplicar a definição. Temos s(0, 3) = fè Iy(tldt = fo? V2dt = E 
Observação 4.7 Se interpretarmos o caminho q(t) como o movimento de uma partícula 
no plano, então a norma (ou módulo) do vetor velocidade |y (t)| é conhecida como veloci- 
dade escalar. E assim, com esta interpretação, o comprimento de arco é a distância total 
percorrida pela partícula ao longo do caminho q. 


Exemplo 4.8 Encontre o comprimento de arco do caminho r(t) = (t2,t3,0) entre t = 0 
et=2. 
Solução: r'(t) = (2t, 3t?, 0), portanto |r'(t)| = V4t2 + 9t4 e 


2 2 
s(0, 2) = f Vv 4t? + 9t4 dt = / V 4 + 98% dt 
0 0 


Fazendo a substituição u = 4 + 9t?, temos du = 18 t dt. Além disto t = 0 corresponde a 
u = 4 e t = 2 corresponde a u = 40, portanto, 


40 du 
0,2) = —. 
Integrando, obtemos 
du PANE, 
[vee a 
Ou seja 
(0,2) = H (403 — 43) = Èpo? — 1] 
Soe a 27 
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Observe que Beate da definição de comprimento de arco e do Teorema Fundamental do 
Cálculo que 4% = |y'(t)|. Novamente, interpretando o oe como a descrição do movi- 
mento de uma partícula percorrendo uma curva C, então E = |y (t)| é a velocidade escalar 
da partícula. Portanto |y'(t)|dt é o produto da ecd. escalar pelo tempo. 

Isto motiva a seguinte definição: 


Definição 4.2.5 Chama-se elemento de comprimento de arco de um caminho y à ex- 
pressão ds = |y (t)|dt 


4.3 Integral Curvilínea de uma função escalar 


Definição 4.3.1 Seja f : U C R? > R uma função contínua e y : [a,b] > R3, y(t) = 
(x(t), y(t), (z(t)) com imagem C, então a integral curvilínea da função escalar f em C é 


definida como: 
[ fds = ia f(y (t)|dt 


Se y(t) = (x(t), y(t), z(t) então |y (t)| = yx (t) 2 + 2!(t)? e escrevemos 


Jis- f reo t)) z (t) )2 + z'(t)2 dt 


Exemplo 4.9 Encontre o valor da integral fẹ yds para C uma curva imagem de y(t) = 
Ci) ts, 
Solução: y'(t) = (¢, 2, —2t) , | (E| = Vt? + 4 + 4t?, y(t) = t, portanto, 


1 
[tas [ t V52+4 dt 
C 


0 


Fazendo a substituição de variáveis u = 4 + 5t? temos du = 10t , quando t = 0, u = 4 
quando t = 1, u = 9. Logo 


1 1 9 ï 3/29 2 
[tas [ne ria == [ pe a 
E 1 


A 5 3h 15 


Um caso particular importante da Integral Curvilinea de uma função escalar ocorre quando 
C é uma curva plana. Vamos examinar este caso com detalhes. Suponha que y(t) é um 
caminho tal que a imagem C é uma curva plana cujos pontos estão no plano xy. Seja f 
uma função real de duas variáveis. A Integral Curvilínea da função escalar f é dada então 
por 


[ten is= f sel, t)) Jel)? + VE de. 


Quando f(x,y) > 0, esta integral possui a interpretação geométrica da “área de uma 
cerca”. Para ver isto imagine a imagem C da função y(t) como como base da cerca e 
para cada (x,y) € C imagine f(x,y) como a altura da cerca no ponto (x, y)(veja a figura). 
Vejamos um exemplo desta situação. 
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f(x), yO) 


Exemplo 4.10 Deseja-se construir uma peça metálica que tem a forma da superfície do 
cilindro x? + y? = 4, compreeendida entre os planos z = 0 e z = 4— 2 — y. Se o metro 
quadrado da peca custa R$ 25,00, calcule o preco total da peca. 

Solução: A base da peça é dada pela circunferência parametrizada por y(t) = (x(t), y(t)) = 
(2cost,2sent),0 < t < 27 e a altura da peca metálica em cada ponto (x,y) € C ima- 
gem de y(t) é dada pela função f(x,y) = 4— x —y. y (t) = (—2sent,2 cost) e portanto 
Iv (£)|| = V4 = 2. Temos então 


27 
OLMO! a= [ (A Dees ~Saen de 1. 
C 0 


O Custo total será portanto 25 x 16a = 4007 reais. 


zh 


y 


Veremos outra aplicações da integral curvilínea de uma função escalar e também um outro 
conceito importante a Integral Curvilinea de um campo vetorial. Mas antes a faremos 


algumas observações uma pouco mais teóricas que nos ajudam a entender o conceito de 
integral curvilínea. 


Observação 4.11 A definição de integral sobre uma curva é motivada do seguinte modo: 
para definir a integral de uma função contínua f em uma curva C que é a imagem de 
um caminho regular, y : [a,b] — R a ideia é seguir os passos da definição da integral de 
Riemann em uma variável. 

Primeiramente, tomamos uma partição Q = UC; da curva em pequenos pedaços e calcu- 
lamos o valor de f(p;) em um ponto p; € C; em cada um desses pedaços. Em seguida, 
toma-se a soma de s(f,Q) = >> f(p;)L(Ci;), para L(C;) o comprimento do pedaço Cj. 
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Observe que, se a correspondência t ++ y(t) é injetiva, então uma partição O corresponde 
a uma partição no dominio a =to<ti<..<tr=beY(t;) = pi. E o comprimento de 
cada pedaço L(C;) é aproximado por |y'(t;)|A;t. De modo que escrevemos 


D=5 J)O) => fat) (last. 


Finalmente, tomamos o limite limjg)_,9 s(f, Q). Observe que 37 f(y(t;))|y (t;)|Ast é pre- 
cisamente uma soma de Riemann para a função f(y(t))|y (t)|. De modo que o limite 


lim¡o¡>0 s(f, Q), se existir, é igual a p FAOI ldt. 


Observação 4.12 A integral de uma função sobre uma curva não depende da parame- 
trização da curva: suponha que t € [a,b] > y(t) e u € [c,d] + ó(u) sejam dois caminhos 
regulares que possuem a mesma imagem, C com as mesmas extremidades, y(a) = ó(c) e 
y(b) = 6(d). 

Se existe uma mudança de parâmetro t +> u(t) tal que q(t) = d(u(t)). Então, pela Regra 
da Cadeia, y(t) = 6’(u(t))u’(t) e |y E| = [8 ult) lu (t), pois u'(t) 4 0 ( por que?). 


b 
Pucoroja= f EONDU Oat 


Mas aplicando a fórmula da mudança de variáveis em integral de uma variável, temos 
du = u'(t)dt, a = u(c) e b = u(d) e 


f EOD O= MCO CO 


Concluímos assim que a integral de uma função sobre uma curva não depende do caminho 
(ou parametrização) que a descreve. Este é o sentido da expressão integral curvilínea: 


fas 


Observação 4.13 Se y(t) é um caminho regular, então existe uma função T(s) tal que o 
caminho r(s) = y(T(s)) tem vetor tangente de norma igual a 1, isto é: |r(s)| = 1. Neste 
caso, dizemos que o caminho está parametrizado pelo comprimento de arco. 


Exemplo 4.14 Verifique se o caminho y(t) = (2 cos(t), 2sen (t)) está parametrizado pelo 
comprimento de arco. Em caso negativo tente achar uma mudança de parâmetro de modo 
que o novo caminho esteja parametrizado pelo comprimento de arco: 

Solução: Iniciamos com o cálculo do vetor velocidade: 


y (t) = (—2sen (t), 2 cos(t)). 


Portanto, |y (t)| = 2 e y nao está parametrizado pelo comprimento de arco. Entretanto, 
se definimos um novo parámetro t = 5, o caminho a(s) = (2cos(5), 2sen (5)) tem vetor 
velocidade igual a a(s) = (—sen (5), cos(5)) que satisfaz |a’(s)| = 1. Ou seja, a está 
parametrizada pelo comprimento de arco. 
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Vejamos como obter teoricamente a função T'(s) no caso geral: 
Primeiramente observe que se fixarmos um ponto tọ no domínio do caminho y então, o 
comprimeto de arco de Y(to) até y(t) é uma função de t definida por: 


T 


(2) = s(to,2) = [ Old (4.1) 


to 


O Teorema Fundamental do Cálculo nos diz que 4 = |7(t)|. 

Como y é um caminho regular, então |y'(t)| > 0. 

Ou seja, a função que a cada T associa o comprimento de arco s(T) possui derivada 
positiva. Pode-se concluir que s(T) é uma função crescente de T ( por que?). 

Usando a interpretação do caminho como a descrição do movimento de uma partícula no 
espaço, isto deveria ser claro, pois, se a velocidade escalar é positiva então, quanto maior 
o tempo de percurso, maior será o comprimento de arco percorrido no caminho. 

Sendo s(T) uma função crescente, então possui uma inversa, que denotamos por T(s). 
Isto significa que (teoricamente) podemos descrever os pontos do caminho usando o com- 
primento de arco como parâmetro. Ou seja um ponto é localizado pelo comprimento de 
arco percorrido a partir de um ponto fixado. 

Além disso, o Teorema da Função Inversa nos garante que a função T(s) é diferenciável 
com 


aT 11 
ds $ WT) 


Se r(s) = y(T(s)), entáo usando a regra da cadeia, 


(= ENE = VE ae 


Portanto, 


Como queríamos demonstrar. 


4.4 Campo de vetores 


A idéia de campos de vetores vem da Física e corresponde a associar a cada ponto do 
plano ou do espaço uma grandeza que possui direção, sentido e tamanho (módulo), ou 
seja, um vetor. Por exemplo, a velocidade de uma partícula movendo-se no espaco é uma 
grandeza vetorial, pois possui direção, sentido e módulo (ou velocidade escalar). 
Sabemos que os vetores do plano e = (1,0) e eg = (0,1) formam uma base de R?. 
Recordando o que isto quer dizer: cada vetor do plano se escreve de maneira única como 
combinação linear de e, e e». Em símbolos, Vu € R? existe um único par de números reais 
a e b tais que v = ae; + bes. 


Definição 4.4.1 Um Campo de Vetores no plano é uma aplicação que a cada ponto (x,y) 
do plano faz corresponder um vetor X(x,y) = F(x,y)e1 + G(x,y)ez ou simplesmente 
X(x,y) = ((F(z, y), G(x,y)) quando a base estiver escolhida. Utilizaremos também a 
notação X quando quisermos enfatizar que para cada ponto (x,y) do plano o campo X 
associa um vetor. 
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Um campo de vetores é dito diferenciável (respectivamente de classe C*) quando cada uma 
das funções F(x,y) e G(x,y) for diferenciável (respectivamente, tiver todas as derivadas 
parciais contínuas). Nestas notas, exceto quando se supõe o contrário, iremos trabalhar 
com campos de classe C?. 


Exemplo 4.15 Damos a seguir 3 exemplos de campos de vetores no plano: 
1. X(x, y) = (2x, 3y). 
2. X(x, y) = (zy, x? — y’). 
3. X(x, y) = (ae Bag) com domínio o conjunto R? \ {(0,0)}. 


Exemplo 4.16 Campo gradiente: se U(x,y) é uma função diferenciável então o campo 
gradiente de U é definido por X(x,y) = VU(x,y) = (&, ow). O campo gradiente é 
bastante estudado no Cálculo Diferencial de Várias Variáveis e uma de suas propriedades 
mais importantes é que aponta sempre na direção de crescimento máximo da função (Ver 
também a Seção 4.7). 


Exemplo 4.17 Igualmente importante é o Campo Hamiltoniano definido por Xy (x,y) = 
OH OH 

(Sy ar) 

Os Campos Hamiltonianos são estudados na Mecânica Clássica. Por exemplo, se H(z, y) = 

my? +U (x), é do tipo energia cinética mais energia potencial, então o Campo Hamiltoniano 


associado se escreve Xy (z, y) = (my, —U'(x)). 


A noção de campo vetorial estende-se naturalmente ao espaço tri-dimensional: é uma 
aplicação que a cada ponto (x,y,z) € Rê associa um vetor 


Xx, y, 2) = (F(x,y, 2), G(x, y, 2), A (a, y, 2)). 


Alternativamente, usando a base canónica de R? dada pelos vetores e, = (1,0,0), e2 = 
(0, 1,0) e ez = (0, 0,1), todo campo vetorial se escreve X (x,y, z) = F(z, y, Je +G(a, y, 2)e2+ 
H(x,y,2)es. Também é usada comumente a notação i, j e k para os vetores da base 
canónica de modo que um campo de vetores se escreve X (x, y, z) = F(x, y, z)i+G(a, y, z)j+ 
H(x,y,z)k. Utilizaremos também a notação X quando quisermos enfatizar que para cada 


ponto (x,y,z) do espaço o Campo X associa um vetor. 


Exemplo 4.18 Campo gravitacional: segundo a Lei Gravitacional de Newton, a forca 
gravitacional exercida por um objeto de massa M colocado na origem, sobre um outro 
objeto, de massa m, situado na posição (x,y,z) é igual a 


mMG 
X(z,y, z) = r3 (x,y,z) 


onde G é a constante de gravitação universal e r = 12 + y? + 22 é a distância do objeto 
até a origem. Observe que o campo X não está definido na origem (por que?). Uma 
expressão similar é dada para um campo elétrico que é dado pela força exercida por uma 
carga pontual Q na origem sobre outra q na posição (x,y, z): 


eqQ 
Y (x.y, z) = “3 (0, Y z). 
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4.5 Integral Curvilinea de um campo vetorial 


Se X é um campo vetorial no espaço então uma partícula neste espaço (por exemplo uma 
massa em um campo gravitacional) vai ser submetida à força X. Suponha que a partícula 
se movimenta ao longo de uma curva C sob a ação de uma força X. Um dos conceitos 
fundamentais da física é o trabalho realizado por X ao longo de C. Veremos que este 
trabalho é medido por uma integral sobre a curva. 

Inicialmente, suponha que a trajetória da partícula é um vetor deslocamento AB eo 
campo é constante. Neste caso sabemos que o trabalho é dado pelo produto escalar do 
campo pelo vetor A 


X. AB= (força ) x (deslocamento na direção da força) 


De uma maneira mais geral, se o caminho é um caminho curvo no espaço, podemos 
imaginar que ele é constituído por uma sucessão de deslocamentos retilíneos infinitesimais. 
Da mesma forma que procedemos na dedução da fórmula para a integral curvilínea de 
uma função escalar, ver a Observação 4.11, chegamos a seguinte fórmula para o cálculo do 
trabalho realizado por um campo X(x,y,z) no espaço sobre uma partícula que percorre 
um caminho regular y : [a,b] > R3, y(t) = (x(t), y(t), z(t): 


b 
trabalho realizado por X = l X (y(t)).9/(t) dt 


Não vamos fazer uma dedução rigorosa desta fórmula mas podemos justificá-la da se- 
guinte maneira. Primeiramente, tomamos uma partição Q do intervalo [a,b], dada por 
a = to < ...ti < tn = b e obtemos uma poligonal no espaço com vértices y(ti) = 
(x(t), y(ti), 2(t)),0 < i < n(ver figura). 


/ 


y DAL 


Se t varia em um intervalo de t; para t;+1 que vamos supor pequeno o deslocamento da 
partícula é aproximado pelo vetor As; = y(t;+1) — y(t;) e o campo X pode ser pensado 
como constante e igual a X(y(t;)) no intervalo [t;,t;+1]. Supondo que y'(t) existe para 
todo t € [a,b] então pela definição de derivada, temos que 


As; & y (t)At; 


Concluimos que o trabalho realizado para deslocar uma partícula de y(t;) para y(ti+1) é 
aproximadamente 


X (y(ti) Asi = X(y(t;)).y (t;)At; 


AULA 4: INTEGRAL CURVILINEA 


63 


64 


Assim o trabalho W realizado pela força X para deslocar uma partícula ao longo de C é 
aproximadamente: 


s(X,Q) = Y X(y(t))As; => XY) (a) Ait. 


Finalmente, tomamos o limite limjg)-,9 s(X, Q).Observe que quando |Q| => 0 a apro- 
ximacao se torna cada vez melhor e que portanto é razoável tomar como nossa definição 
de trabalho a integral 


b 
w= xr 
a 
Esta nocáo de trabalho, fundamental na física, nos leva a definir a integral curvilínea de 


um campo vetorial da seguinte maneira. 


Definição 4.5.1 Seja X um campo vetorial no R? contínuo no caminho regular y : 
[a,b] — R3. A Integral curvilínea do campo X ao longo de y é definida da seguinte 
maneira: 


b 
[xo = J X(y(t).y (t)dt 
y a 
Utilizando coordenadas, se o campo vetorial se escreve 


X (x,y,z) = (F(x,y, 2), G(x, y, 2), H (x,y, z)) e yt) = (z), yt), 2), 


X (y(t). (6) = F(a(t), y(t), 2) a") + Glelt), y(t), z(t))y (6) + H (x(t), y(t), 2(t)) 2") 


b 
| xas= [EGO + coro + Homo. 
Y a 


Exemplo 4.19 Considere o campo vetorial X(x,y,z) = (x,y,z?) e o caminho y(t) = 
(cos(t), sen (t),t) onde 0 < t < 7. Calcule a integral curvilinea i X.ds. 

Solução: O campo ao longo do caminho é igual a X(y(t)) = (cos(t), sen (t), t?) e o vetor 
tangente ao caminho é igual a y'(t) = (—sen (t), cos(t), 1) para 0< t< 7. 

Logo, X(y(t)).y'(t) = t? e, por definição, 


g T 


T t3 
[xa= | t dt = — 
y 0 3 10 


Como já vimos, o conceito de integral curvilínea de um campo está associado ao conceito 
de Trabalho realizado ao movermos uma partícula ao longo de uma curva a, sob a influência 
de um campo de forças X = (F,G). Considere o seguinte exemplo: 


Exemplo 4.20 Encontre o trabalho realizado pelo campo X (x,y) = (xy, y?) ao se mover 
uma partícula ao longo de uma curva C descrita pelo caminho r(t) = (3 cos(t), 2sen (t)) 
para0<t<r. 

Solução: X(r(t)) = (6cos(t)sen (t), 4sen?(t)) e r’(t) = (—3sen (t), 2 cos(t)). 

Logo X(r(t)).r'(t) = —18 cos(t)sen 2(t) + 8 cos(t)sen 2(t) = — 10 cos(t)sen ?(t). 

Portanto Wc = Ju X.dr = Jo [-10 cos(t)sen 2(t)]dt = cost (8), = 5 
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Dbservação 4.21 A componente de X na qe do vetor unitário da tangente a y no 
onto t é o produto escalar f(t) = X(y(t)). rue ae : que é uma funcáo escalar do parámetro t. 
[sto quer dizer que ao projetarmos ortogonalmente X sobre o vetor unitário da tangente 
bteremos a função escalar f(t). Obtemos assim uma relação entre a integral curvilínea 
lo campo e a integral de f no caminho y que pode ser muito útil: 


[soa = f xow O oja= xa) (tdt 


Observação 4.22 Seja y(t) = (x(t), y(t), z(t)) um caminho regular. Suponha que g seja 
ima função diferenciável com g'(u) > 0 tal que g(ua) = a,glup) = b. Definimos um 
aovo caminho regular a(u) = y(g(u)) obtido do caminho anterior por uma mudança de 
parâmetros t = g(u). Pela Regra da Cadeia a’(u) = y (g(u))g'(u). Portanto, 


Por definição, dado um campo X 


f k= / o Aaa 
[ xas= [xy [IMA wan 


Aplicando a fórmula de mudança de variáveis nesta útima integral, obtemos: 


Up b 
[xa= | Xy [xt rs | Xa. 
a Ua a y 


Em outras palavras, a integral curvilinea de um campo X ao longo de um caminho, não 
depende da parametrização do caminho. 
Faz sentido, portanto, definir a Integral Curvilínea de X sobre a curva regular C, o traço 
‘ou imagem) de um caminho y regular: 


Logo, 


/ Pe Lea ie E X(-(t)).7 (tat 
C Y 


A expressão Fdz + Gdy + Hdz deve ser interpretada como o produto escalar do campo 
vetorial X = (F,G,H) com o vetor tangente dado por uma parametrizacáo regular y da 
curva C. A integral de linha na forma 
1 X.dr. 
C 


$ X.dr 
C 


dara a integral na curva completa (uma volta completa). 


Quando a curva é fechada escreve-se 
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Observação 4.23 Se y(t) = (x(t), y(t), z(t)) e B(t) = y(a+b—t), então o traço de ambos 
os caminhos coincide, mas 3(a) = y(b) e 6(b) = y(a). 
Isto é os sentidos dos percursos são contrários: B'(t) = —y (a +b-—+t). 


b b 
ds = a a (t))dt = a+b . (a+b d 
J x I X(a(t)).a!(t))dt f EE t))dt 


Fazendo a mudança de variáveis, u = a + b — t, du = —dt, 


b a 
[x=] Xola+b- 0a +o- e)at= | ss 


- xo) yan = - J xao 


Ou seja, ao mudarmos o sentido do percurso de uma curva C, a integral do campo vetorial 


muda de sinal: 
f Xas=- | Xas 
B Y 


Isto significa que se denotarmos por —C' o traço de y percorrido no sentido contrário ao 
de C, então 


f Fdz+Gdy+Hdz=- | Fdz + Gdy + Has, 
=C C 


Exercício 4.24 Verifique esta observação no seguinte caso: X(x, y, z) = 2yi+ xzj + yzk, 
C a imagem de q(t) = (t, t?, t3) e —C a imagem de y(—t) = (—t,t2, —tº), para —1 < t < 2. 


Exemplo 4.25 Geralmente náo se faz mencáo explícita ao campo de vetores e escrevemos 
uma integral curvilínea na forma de Fdx + Gdy + Hdz. Se C é o traço do caminho 


y(t) = (t,0,0) , a < t < b, então fo Fdz + Gdy + Hdz = J? F(t,0, O)dt. 


Exemplo 4.26 Calcular fo 2ydx + zdy + vzdz, C = (a? + y? = 1,z = 0} parametrizado 
no sentido anti-horário ( positivo). 
Solução: observe que y(t) = (cos(t), sen (t),0) é uma parametrização de C, 0 < t < 27. 


27 
/ 2ydx + xdy + xzdz = / [2sen (t)(—sen (t)) + cos(t) cos(t)]dt. 
C 0 


Usando identidades trigonométricas, cos(2t) = cos*(t) — sen?(t) e cos?(t) + sen?(t) = 1 
temos sen (t)(—sen (t)) + cos(t) cos(t) = 2cos(2t) — 1 e substituindo na última integral 
obtemos: 


27 


27 
/ ydz + zdy + vzdz = / [2 cos(2t) — 1]dt = [—sen (2t) — t] a —27. 
C 0 


Pergunta: qual seria o resultado se a curva estivesse orientada no sentido contrário? 
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Exercício 4.27 Sejam y(t) = (e~’sen(t),e~'cos(t),t) para O < t < 3 um caminho e 


X(xz,y,z) = (EG H) = sopro (ty, z). Encontre a integral curvilinea: 
[ras + Gay + Hdz 
fi 


Definição 4.5.2 A integral curvilinea pode ser estendida para um caminho regular por 
partes, ou seja um caminho continuo  : [a,b] > R? constituído de uma união finita de 
caminhos regulares que se intersectam em no máximo um ponto. 

Portanto existe um número finito de pontos a = ay < a, < ... < ak = b tal que a restrição 
de y a cada um dos sub-intervalos [a;, a;+1] é um caminho regular. 
Se denotamos por C o trago de y e por C; o traço do caminho Ylļja;, 
sub-intervalo [a;, a;+1], então, definimos 


a;+1)» restrição de y ao 


k 
L Fdx + Gdy + Hdz = >| Fdz + Gdy + Hdz 
c i=1 64 


Exemplo 4.28 Calcular $, zdx + zdy + ydz para C o triângulo de vértices O = (0,0,0), 
A = (1,1,0) e B = (1,1,1), percorrido na seguinte ordem: de O para A, de A para B e 
de B para O. 

Solução: C é um triângulo formado por três segmentos (arestas): 

Cı, de O para A: y1 (t) = (t,t,0), para 0< ¢ < 1; 

Ca, de A para B: y2(t) = (1,1,t), para0<t<le 

C3, de B para O: y3(t) = (1—t,1-—#,1-1t), para 0 < t < 1. 

Observe que estamos usando o fato de que a integral curvilinea nao depende da parame- 
trização. 

Logo fo zda + zdy + ydz = o Jo, 2dz + zdy + ydz. 

Calculamos separadamente 


1 
Jc, “dz + ady + ydz = ht + t)dt = el =1 


g 1 
Jo, z4x xdy + ydz de dt Ho 1 


Jo, zdx + zdy + ydz = J —3(1—t)dt = 3 (1 — t)? 


Portanto, 


3 1 
$ ade + ady +ydz=1+1- 5 = 
C 2 2 


Exercício 4.29 Encontre a integral curvilínea bo aydx + xy*dy para C o triángulo de 
vértices A = (0,0), B = (0,1) e C = (1,1), orientado no sentido: de A para B, de B para 
C e de C para A. vskip 5cm 


4.6 Teorema de Green 


Nesta seção estaremos restritos a caminhos e campos vetoriais no plano. O objetivo 
desta seção é expor o Teorema de Green, que estabelece uma relação entre a integral 
curvilínea em um caminho fechado e a integral dupla de uma função ( o ”rotacional” do 
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campo) na região limitada pelo traço do caminho. Veremos que de certa forma estaremos 
generalizando o Teorema Fundamental do Cálculo. 

Seja R uma região plana limitada por uma curva y união de um número finito de curvas 
disjuntas y;, para i = 1,2,...,n. Cada uma das curvas y; é regular por partes e está 
orientada de modo que ao percorrê-la, o interior da região está sempre à esquerda. Para 
enunciar o Teorema de Green vamos definir o conceito do rotacional de um campo vetorial 
plano. 


Definição 4.6.1 Dado um campo vetorial X no plano, o rotacional de X, é a diferença 
entre a derivada parcial da segunda componente do campo em relação à primeira variável 
e a derivada parcial da primeira componente em relação à segunda variável: 


rot(X) = [Ge(#,y) = Fyle= === 50" 


Teorema 4.6.1 (Green) Seja X(x,y) = (F(x,y),G(x,y)) um campo vetorial definido 
em um subconjunto aberto U do plano tal que RC U, X : U —> R?, onde R é uma região 
plana limitada por uma curva orientada y como acima. 

Então 


$ Fla yde + G(e.a)dy = | [1Gs(0.9) — Fyle y)]dA 


Em palavras, para campos vetorias X no plano, O Teorema de Green afirma que a integral 
curvilínea em uma união disjunta de curvas regulares por partes que limita uma região 
plana R é igual à integral dupla do rotacional do campo em R. 

Antes de fazer a demonstração damos a seguinte definição: 


Definição 4.6.2 Um região R do plano é chamada uma região simples se ela pode ser 
descrita simultaneamente com uma região do Tipo I e como uma região do tipo II (ver 
Seção 1.4). 


Demonstração: Suponha inicialmente o caso particular em que R seja uma região sim- 
ples. Ou seja, vamos supor que R = {a < x < b3gi(x) < y < go(x)} e simultâneamente 
R= {c < y < d; hi (y) < £ < ga(y)). 


A 
y 


ye 


LIÇÕES DE CÁLCULO INTEGRAL EM VÁRIAS VARIÁVEIS 


Desta forma, a curva y que limita R (ou fronteira de R) é descrita de duas maneiras: A 
primeira: y = y U y3 com y(x) = (2, gı(x)) e y(x) = (a +b — x,g(a+b-— x)) para 
a <x <b. Observe que as curvas e y; e y3 são percorridas em sentidos contrários. 

A segunda: y = y2 U ya com y2(y) = (gı(y) y) e yaly) = (ha(c + d — y,c + d — y) para 
c<y<d. A curva ya é percorrida no sentido positivo do eixo Oy, enquanto a curva yq é 
percorrida em sentido contrário. 

A prova do teorema é consequência do Teorema Fundamental do Cálculo. Calculemos 


| [Cen -Fesa | f Geaa- f f Etemia 


Para a primeira integral, escolhemos integrar primeiramente em relação a x e usamos a 
seguinte descrição para R : R = {c < y < d; hı (y) < xz < go(y)}. 
Assim, pelo Teorema de Fubini, 


| f Galena = Flo G, (z, y)dajdy. 


Pelo Teorema Fundamental do Cálculo: 


d 
J | Gelenaa= feat). — Gn) wld. 
R c 


Mas o lado direito da equação é igual a 


Geva- f G(x, y)dy. 
Ya 


y2 


JJ, esmas f cmd | Ets 


Para a segunda parcela, integramos primeiramente em relacáo a y: 


| ft ias [fe r(x, y)dylde. 
| | Frenda = f Eeo) Gatas 


Mas o termo direito da equação é igual a 


EA Plajdo = f Cds 


J |Een f Fley)de+ f Pleno 


Subtraindo as igualdades encontramos: 


| fre (£, y) — F(z, y)JdA = 


Portanto, 


Ou 


Substituindo 
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=f cmd | Steady | Fæy)de+ | F(x, y)dz. 


13 
Concluimos assim que 


f fico) — F,(x,y)|dA = [Feudo + Ga, y)dy. 


Y 
Para o caso geral, basta decompor a região R em sub-regiões simples R; ; como acima por 
meio da introdução de segmentos verticais ou horizontais. 


Aplica-se o resultado em cada sub-região R;; e, em seguida, observa-se cada segmento 
introduzido é percorrido duas vezes em sentido contrário de modo que ao somarmos 


Df, Gee — Ele, yJdA 


as integrais curvilíneas correspondentes se cancelam, restando apenas a integral curvilínes 
sobre os arcos de curvas ;,2 = 1,2,...,n da fronteira de R, concluindo assim a demons- 
tração do Teorema de Green. 


Exemplo 4.30 Seja y uma curva regular simples e fechada limitando uma região R. Seja 
X(x,y) = (x,y) um campo.Entáo 


fede +yay = f [0 dA = 0. 
Y R 


Exemplo 4.31 Seja novamente y uma curva regular simples e fechada limitando uma 
região R. Considere agora o campo X(x,y) = (—y, x), então: 


$ -vaz + zdy = T 2dA = 2 área (R). 
y R 


Exercício 4.32 Seja y uma curva regular simples e fechada limitando uma região R. Se 
Fy, (x,y) = Gala, y), V(x,y) € R, calcule 


f F(a,y)dx + G(x, y)dy 
J 


O Teorema de Green tem muitas aplicações na Matemática. Entre as mais importantes 
que estudamos no curso de graduação citamos: a integral de funções complexas e na 
Geometria Diferencial. Também encontramos aplicações do Teorema de Green na Fisica 
quando estudamos os chamados campos conservativos. É o que veremos na próxima seção. 
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4.7 Campos Conservativos no Plano 


Seja r(t), a < t < b uma parametrização regular por partes de uma curva orientada C C R° 
e X um campo continuo sobre C. Então vimos a seguinte definição: 


[ xare f xeo 


Em geral a integral curvilinea ao longo de uma curva C parametrizada por r(t),a <t<b 

depende da curva C e dos pontos inicial e final da curva A = r(a) e B = r(B). Veremos 

agora um caso particular importante em que a integral curvilínea depende apenas dos 

pontos inicial e final A e B mas não da curva C. Para isto, considere o campo X = Vf, o 
mer a — (of af Of 

gradiente da função f, dado por Vf = (55, By) Aa) 


Teorema 4.33 Seja X um campo vetorial tal que exista uma função f satisfazendo 
Vf = X. Se C C RÌ é uma curva regular por partes com pontos inicial e final A e B 
respectivamente, então temos: 


fxar= | Sidr = iB- A 


Demonstração: Como r(a) = A e r(b) = B temos 


b 
/ ALS = f Vif (r(t)).r'(t)dt 
C a 
Seja g(t) = f(r(t)) a<t<b. 
Temos que 
(= Ofdx 0fdy | Of dz 
4 dx dt Oydt ` Əz dt 


pela regra da cadeia. Mas então segue do Teorema Fundamental do Cálculo que 


= Vf (r(t))1(0), 


b 
[ Xar= [ g(t) = 9(b) — g(a) = F(r(b)) — F(r(a)) = f(B) — HA). 


Definição 4.7.1 Um campo X que satisfaça o teorema anterior, isto é tal que exista uma 
função f com X = Vf é chamado um campo gradiente. 


O seguinte teorema caracteriza os campos gradientes no plano. 


Teorema 4.7.1 Seja X = (F,G) um campo vetorial no plano com componentes contínuas 
e diferenciáveis em uma aberto U C R?. Então são equivalentes: 


1. bo X.dr = 0 para toda curva fechada C regular por partes contida em U. 


2. A integral curvilínea de X do ponto A até o ponto B independe da curva que liga A 
a B. 
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3. X é um campo gradiente de alguma função f em U. 


4 OG _ OF 
Ox — Oy’ 
Demonstração: 


Vamos mostrar que (1) => (2) => 3 => (4) => (1). Inicialmente para ver que (1) =; 
(2), considere duas curvas C4 e C2 regulares por partes contidas em U e unindo os pontos 
A e B (ver figura). 


A 


Pela Observação 4.23, dada a curva C2, existe a curva —C uma curva idéntica a C2 mas 
com ponto inicial B e ponto final A. Então C = C1 U —C5 é uma curva fechada e por (1) 
temos que: 


$ X.dr = X.dr — X.dr = 0. 
C Ci C2 


Segue que e, X.dr = Jos X.dr. 

Em seguida vamos mostrar que (2) => (3). Seja C uma curva ligando o ponto (0,0) a 
um ponto (x,y) e suponha que C possui uma parametrização r(t). Defina f = Ja X.dr. 
Por hipótese f independe de C. Vamos mostrar que f = VX. Para isto escolha o caminho 
C = C1 U Co onde Cı é parametrizada por rı = (t,0), 0<t< x e Ca é parametrizada 
por rg = (x,t), 0<t< y. Temos: 


ale [Proa 4 | cena 


Segue que of = G(x,y). De maneira semelhante permutando x e y obtemos que = = 
F(x,y), concluindo esta parte da demonstração. 


Para ver que (3) => (4), suponha que X = Vf, para alguma f. Isto significa que F = of 


2 2 
eG= a Mas então = = E e e = ZE. Mas sabemos do curso de cálculo diferencial 


2 
em várias variáveis que EE = = a nestas condições, concluindo a demonstração desta 


parte. 


Finalmente é fácil mostrar que (4) => (1), pois o resultado segue diretamente do Teorema 


de Green, pois 
pxar= [| (cg y — Fr) ) dedy= ff 0 dx dy = 0, 


onde D é o interior do caminho fechado C. 
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Definição 4.7.2 Uma função f que satisfaz as condições equivalentes do Teorema 4.7.1 
é chamada uma Função potencial, ou seja, f é potencial se existe um campo X tal que 
f =VX. Neste caso dizemos que o campo X é um Campo Conservativo 


Observação 4.34 Podemos usar o Teorema 4.7.1 para encontrar uma função potencial de 
um campo conservativo e para calcular facilmente algumas integrais curvilíneas de campos 
conservativos. Observe o exemplo a seguir. 


Exemplo 4.35 Considere o campo X = (F,G) = (2xy, 1? — y?), no plano. 


1. Verifique que X é um campo conservativo. 


Solução: Como as condições do Teorema 4.7.1 são equivalentes basta mostrar que 
a condição (4) é satisfeita. Com efeito Fy = 2x e Gy = 2x. Logo Gr — Fy=0e0 
campo é conservativo. 


2. Encontre uma função potencial para o campo X. 


Solução:Estamos procurando uma função f tal que Vf = X = (2xy,x? — y?). 
Sabemos que of = 2xy. Integrando em relação a x podemos concluir que f(x,y) = 


ay + hi(y) onde h não depende de x. De maneira semelhante, temos que i = 


r? — y?, donde concluimos que f(x,y) = x2y — e + ho(x), onde h(x) não depende 
de y. Comparando as expressoes obtidas para f encontramos uma funcáo potencial 
f = ry- r, Observe que é esta é apenas uma das possibilidades, não existe 
unicidade na função potencial procurada. 


3. Calcule Jo X.dr onde C é o arco de parábola dado por (t,t?) para 0 <t <1. 


Solugao:Como encontramos a função potencial do campo temos 


f xar= 108) E 
onde A = (0,0) e B = (1,1). 

Exercício 4.36 Encontre outras funções potenciais para o mesmo campo X. 
Observação 4.37 Se rot(X) = Gr — Fy = 0 então o trabalho realizado não depende 


da curva, apenas do ponto incial e do ponto final. Isto é uma consequência imediata do 
teorema acima. 
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4.8 EXERCICIOS 


il. 


10. 


11. 


Verifique se cada um dos caminhos abaixo é regular ou nao: 


(a) y(t) = #0, #4) 

(b) y(t) = Gano. te). 0<t< r 
(c) y(t) = (0,1%, 12) 

(d) y(t) = (etcos(t), e’sen(t), t) 


. Dados dois pontos A = (a1, a2,a3) e B = (b1,b2,b3) em R$, encontre um caminho 


diferenciável cuja imagem é a reta que passa pelos pontos. 


. Se N = (1,0), S = (—1,0), encontre um caminho cuja imagem está contida no 


círculo x? + y? = 1 e se inicia em N e termina em S. 


. Mostre que o caminho 0 ++ (3cos(0), 4sen(0)) é regular e sua imagem está contida 


2 
na elipse de equação E Le > 1 


. Encontre fydx, C := (x, x), para 0 < z < 1. 


. Deseja-se construir uma peça metálica que tem a forma de um cilindro x? + y? = 9, 


compreendida entre os planos z = 0 e z = 9 — x — y. Se o metro quadrado da peça 
custa R$ 50,00 calcule o custo total da peça. Faça um esboço da peça. 


. Seja C a fronteira do quadrado no plano xy de vértices A = (0,0), B = (1,0), 


C = (1,1) e D = (0,1) orientada no sentido anti-horário. Considere o campo do 
plano dado por F(x,y) = (x?, xy). Calcule a integral curvilínea (ou de linha) do 
campo F ao longo da curva orientada C. 


(a) Escolha um caminho y ligando o ponto (1, 1) ao ponto (2, —3) e calcule f y 21ydi+ 
x?dy.(Sugestão: utilize um segmento de reta.) 


(b) A integral calculada depende do caminho? Por que? (Sugestão: Aplique o 
Teorema 4.7.1). 


. Encontre a integral curvilinea to x22ydx + vydy para C a curva fechada formada pelo 


arco de parábola y = x, para 0 < x < 1, percorrida no sentido crescente de x e pelo 
arco de parábola x = y?, percorrida no sentido decrescente de x (de 1 a 0). 


Encontre le + xy?) dz + (x2y + y? + 3x)dy, para C a elipse a + a = 1 orientada 
positivamente ( sentido anti-horário). 


Encontre fo 2xydx+ y1 + ytdy para C o triángulo de vértices A = (0,0), B = (0,1) 
e C = (3,1), orientado no sentido negativo, isto é: de A para B, de B para C e de 
C para A. 
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12, 


13. 


14. 


15. 


16. 


Atenção! Devemos ter cautela ao aplicar o Teorema 4.7.1. Seja X (x,y) = a gage) 


(a) Verifique que rot(X) = 0. 


(b) Se a(t) = (cos(t), sen(t)) e B(t) = (cos(t), —sen(t)), O < t < 7 são curvas que 
ligam (1,0) a (—1, 0), é verdade que W (a) = W(B) ? 


(c) O que falha? 
Utilize o Teorema de Green para encontrar as seguintes integrais curvilineas: 


(a) fo vyda + xy?dy para C o triângulo de vértices A = (0,0), B = (0,1) e C 
(1,1), orientado no sentido: de A para B, de B para C e de C para A. 


(b) ¢,ayda + x?dy para C a curva constituída pela semi-circunferéncia 1? + y? = 
4,2 > 0 e o segmento x = 0, —1 < y < 1, orientada positivamente. 


(c) fole” + 3y?)dx + (cos(,/y) — 2xy)dy para C o quadrado de vértices (—1,—1), 
(1,—1), (1,1) e (—1, 1), orientado positivamente. 


Se C é o segmento de reta ligando o ponto (x1, y1) ao ponto (£2, ya), mostre que 


ady — ydx = x1 y2 — 1241. 

C 

Se os vértices de um polígono de n lados, na ordem anti-horária, são (21, y1), 
(22, Y2),..., (Tn, Yn), mostre que a área do polígono é 


1 
A= z (2192 £2Y1) + (T2Y3 L3Y2) + 7 (Oo = Batin) ar Eo e 


Determine a área do pentágono com vértices (0,0), (2,1), (1,3), (0,2) e (—1, 1). 
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Teorema de Stokes 


AULA 5: TEOREMA DE STOKES 


Objetivos 5.1 Os objetivos desta Aula sao: 


e calcular a integral de uma função escalar e de um campo vetorial sobre uma região 
contida em uma superficie (integral de superficie); 


e estudar as propriedades de campos vetoriais por meio do rotacional e da divergéncia; 
e introduzir o conceito de fluxo de um campo através de uma superficie; 


e relacionar a integral curvilinea de um campo ao longo de uma curva com o fluxo do 
rotacional do campo na região da superfície limitada pela curva. 


e estudar um tipo de campo importante: os campos conservativos; 


Nosso próximo objetivo é generalizar o Teorema de Green para campos de vetores no espaço 
tridimensional. As generalizações que descreveremos têm aplicações importantes tanto 
na matemática quanto no eletromagnetismo, ou na mecânica dos fluidos, por exemplo. 
Conforme o ponto de vista adotado, há duas generalizações para o Teorema de Green. 

A primeira considera a integral curvilínea sobre uma curva C C S que limita uma região 
D = r(R) contida na imagem de superfície parametrizada regular S. Lembre-se que o Te- 
orema de Green (no plano) relaciona a integral curvilínea sobre uma curva parametrizada 
y a uma integral dupla de uma certa expressão do campo no interior de R. A primeira 
generalização que estudaremos neste capítulo, relaciona a integral curvilínea em C com 
uma integral de superfície no interior da região D. 

A segunda generalização, que será tema do próximo capítulo, relaciona uma integral sobre 
uma superfície parametrizada regular S que limita uma região no espaço 2 com uma 
integral tripla no interior da região. Nesta generalização, ao invés de integral curvilinea 
teremos uma integral numa superfície e no lugar de integral dupla teremos integral tripla. 
Neste sentido dizemos que esta é uma generalização na dimensão. 


Iniciemos, com a generalização da Integral Curvilínea. 


5.1 Integral de Superfície de funções escalares 


No Capítulo II, aplicamos a integral dupla para calcular a área de uma superfície para- 
metrizada. 
Recordemos a definição de superfície parametrizada regular. 


Definição 5.1.1 Uma superfície parametrizada regular é uma aplicação r : U C R? > R°, 
r(u,v) = (x(u, v), y(u, v), z(u, v)) tal que os vetores 

rulu, v) = (Lulu, v), Yulu, v), zulu, v)) e ro(u, v) = (Lulu, v), y(u, v), Zulu, v)) 
são linearmente independentes para todo (u,v) € U. 


A imagem r(U) = S é chamada superfície regular. A expressão dS = |ru xr,|dA é chamada 
elemento de área da superfície parametrizada. Seja f : W CRY > R uma função contínua 
cujo domínio contém S a imagem de uma superfície parametrizada r. 

Sejam R C U uma região contida no domínio de r e D = r(R). 
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Definição 5.1.2 A Integral de Superfície de uma função escalar f em D é definida pela 
seguinte expressão: 


J | t08- | | Feto) atuo) Au) x rolda 


Esta definição é análoga à definição de integral em uma curva tratada no capítulo anterior. 
De maneira semelhante àquele caso é possível provar, usando a fórmula de mudança de 
variáveis para a integral dupla, que, se D C S é a imagem de duas regiões por parame- 
trizações regulares então o valor da integral de superfície é o mesmo. 

Observe também que se f = 1 então a integral de superfície é igual à área da região que 
foi tratada nas aplicações da integral dupla (ver Definição 2.2.2). 


Exemplo 5.2 Encontre a massa de uma placa D descrita pela imagem do disco u? +v? < 
1 pela parametrização r(u,v) = (u,v, u? + v?), cuja densidade é dada por p(x,y,2) = 
a Observe que esta é a parametrização de uma região do paraboloide z = x? + 


y?. Primeiramente calculamos os vetores geradores do plano tangente em cada ponto: 
ry(u,v) = (1,0,2u) e ry(u,v) = (0,1,2v)) e o vetor normal: ru x ry = (-2u,-—2v,1) e 
|ru X Py] = V4u? + 4v? + 1. Por definição, a massa total é igual a: 


= 1 e/a ge, 
[as [Ga 4u? + 4v? + 1dA 


Para calcular a integral dupla, usamos coordenadas polares: 


V4u2 + 40241 1 T r [ r 
JdA = — dôdr = 27 —— dr 
NA aes o Jo v8r? +1 o v8r? +1 


Calculando a integral por meio da substituição t = 8r? + 1, dt = 16rdr, obtemos: 


1 9 
r 1 T 9 T 
27 —— dr = 2r | — — dt = “vil = — 
/ 8r2 +1 1 16vt 4 1 2 


Exercício 5.3 Calcule a integral de superfície da função f(x,y,z) = y sobre o triângulo 
de vértices nos pontos (1,0,0), (0,1,0) e (0,0,1). (Observe que os pontos estão contidos 
no plano z +y +z = 1.) 


5.2 Integral de Superfície de campos vetoriais 


Seja X(x,y,2) = (F(x,y, z), G(x,y, z), H(x,y,z)) um campo de vetores em R3. Seja S 
uma superficie parametrizada por r(u,v) onde (u,v) € R C R?. Definimos: 


Definição 5.2.1 A Integral de Superfície do campo X na superficie S é definida pela 


seguinte expressão: 
f | xas = al X(r(u,v)).(ty X Py) du du = 
S R 


= | | X2(u.0).00u,0).2(u, 0). sriid 
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A integral de superfície definida acima de um campo X sobre uma superfície S é denomi- 
nada o fluxo de X através de S. 


Observação 5.4 Segue da definição que o fluxo depende da orientação da superfície, isto 
é, se usamos uma parametrização da superfície com vetor normal no sentido contrário, 
então o fluxo tem sinal oposto. 


Vejamos a seguir um exemplo de como se calcula esta integral e, em seguida, a razão desta 
denominação. 


Exemplo 5.5 Seja S a esfera de raio 1 e X(x,y,z)= xi + jy+ kz um campo. Va- 
mos calcular f fg X.dS. Para isto consideramos a parametrização da esfera, utilizando 
coordenadas esféricas (ver 3.3), dada por: 


r(u,v) = (cos(u)sen (v), sen (u)sen (v), cos(v)) 


com 0 < u < 2m e 0 < v < m. Denominemos D a região do plano uv dada pelas equações 
acima. 


Temos que ru = (—sen (u)sen (v), sen (u) cos(v),0) e ry = (cos(u) cos(v), sen (u) cos(v), —sen (v)). 


Consequentemente ru X ry = (—sen ?(v) cos(u), —sen 2(v)sen (u), —sen (v) cos(v)). 
Em seguida calculamos X (r (u, v)).(£u X rv). 
Como X(r(u,v)) = (cos(u)sen (v), sen (u)sen (v), cos(v)), temos 


X (r(u,v)). (£u X ro) = (5.1) 
(cos(u)sen (v), sen (w)sen (v), cos(v)).(—sen (v))(—sen (v) cos(u), —sen (v)sen (u), cos(v)) = 
— (sen (v)) (sen 2(v) cos? (u) + sen 2(v)sen Zu + cos? (v) = 


= —sen (v) 


Finalmente: 


f [xas= f f -w du w= [" f” sen) dv du = —4r 


j 


Observação 5.6 Observe que se r(u,v) = (x(u,v), y(u, v), z(u,v)) é a parametrização 
uma superfície regular D = r(R) com vetor unitário da normal igual a N (u,v) = fi 
então o produto escalar X(x(u, v), y(u, v), z(u, v)).N(u,v) = f(u,v) é a componente do 
campo de vetores X na direcáo do vetor normal á superfície. Se imaginarmos uma família 
de curvas tangentes ao campo vetorial X, então a função f(u,v) pode ser interpretada 
como uma medida de quanto o campo entra (caso em que f(u,v) > 0 ou sai ( respectiva- 
mente f(u,v) < 0) da superfície. Podemos, portanto, calcular a integral de superfície da 
funcáo escalar f(u,v) sobre D. 


Definimos o Fluxo de X sobre a superfície como 


f | xas= | f| Xeer) du dv =. 


Mas temos as seguintes igualdades: 


J | sas= | | Xelo) 20,0). veux rulda 
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Tu 


Jf 1as= | | Xele) to) tno) EE eu x ay de. 


/ ftas= / Í X (æ(u,v), y(u, 0), 2(u,0)) (tu x rv)dA = / [ X.dS 


Ou seja temos duas expressões equivalentes para o fluxo. Vamos explorar este fato mais 
adiante. 


Exercício 5.7 No Exemplo 5.5 calculamos o fluxo de X na esfera com vetor normal 
apontando para dentro, dado por —sen (v). Calcule o fluxo de X na mesma esfera com o 
vetor normal apontando para fora. 


oe _ 2x 2y 2z 
Exercício 5.8 Calcule o fluxo de X(x,y,2) = (y. Faye ayaa) sobre a es- 


fera centrada na origem e de raio 2 com vetor normal apontando para fora. 


Observação 5.9 Para entender melhor a Observação 5.6 acima considere a seguinte si- 
tuação física. Imagine um cano de água de seção circular S e considere um campo X que 
expressa o escoamento de um líquido em cada ponto. O fluxo de X através de S mede 
a quantidade de água que está escoando através de S. Imaginemos X constante em cada 
ponto. Se X é normal ao cano o fluxo é zero significando que não há vasão de água através 
do cano. Se X está na direção da normal de S a vasão é máxima para fora ou para dentro 
do cano dependendo do sentido de X. Veja a figura. 


S 


N Es 
SITUAÇÃO DE VASÃO ZERO 


Exercício 5.10 Faça o esboço do cano e do fluxo do campo X quando a vasão é máxima 
para dentro do cano. 


5.3 Teorema de Stokes 


Nesta seção será exposta a primeira generalização do Teorema de Green de que falamos 
no início do Capítulo. Recordando a conclusão do Teorema de Green é que: 


[pars cay= f f iG. - rea. 
€ R 


Designamos o integrando Gys — Fy como o rotacional do campo X. Esta expressão 
generaliza-se para dimensáo trés, o rotacional de um campo de vetores é um outro ve- 
tor definido por 


LIÇÕES DE CÁLCULO INTEGRAL EM VÁRIAS VARIÁVEIS 


Definigao 5.3.1 Considere 
X(x, Y, z) = (F(a, Y, z), G(x, Y, 2), H(z, Y, z)) 


um campo de vetores em R3. Chama-se rotacional de um campo X ao campo rot(X) 
definido por: 


rot(X) = (H} — Gz, Fz — Hz, Gz E) 


É claro que se X é um campo vetorial no plano então X(x,y,z) = (F(x,y), G(x, y), 0) 
então 


rot(X) = (0,0, Gz — Fy). 


Há uma notação muito sugestiva para o rotacional, usando a notação do produto vetorial. 
_(0 0 ə z 
Se V = (a5 By? 5z) então, 


i j k 
rot(X)=VxX=|2 5 2 
F ĠG H 


Exercício 5.11 Encontre rot(X) para o campo X(x, y, 2) = (x7, y?, 22) 


Antes de enunciar o Teorema de Stokes, vamos discutir um pouco o conceito de orientação 
de uma curva. 

Vimos no capítulo sobre integral curvilínea, que o valor da integral troca de sinal quando 
mudamos o sentido de percurso da curva (orientação). Para generalizar o Teorema de 
Green, é preciso dar sentido ao conceito de orientação de uma curva contida em uma 
superfície. Observe que escolhida uma parametrização regular r(u,v) para a superfície, 
temos definido um campo vetorial normal á superfície: 


Dada uma curva regular fechada C C S, que limita uma região D simples na superfície, 
podemos escolher uma parametrização y(t) para C. 

Observe que N(y(t)) e y'(t) são vetores ortogonais, de modo que N(y(t)) x y(t) £ 0 é um 
terceiro vetor ortogonal a ambos. 
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Definição 5.3.2 Nas condições descritas acima, dizemos que a parametrização y(t) dita 
positiva na superficie se a matriz cujas linhas são y(t), N(y(t)) x y (t) e N(y(t)) possui 
determinante positivo. 

Equivalentemente, se o produto misto y (t).[N (y(t) x VH] x N(y(t)) > 0. 


Uma interpretação pictórica desta definição: a parametrização é dita positiva se ao cami- 
nharmos ao longo de y a região D fica à nossa esquerda. 

Este conceito de orientação coincide com o formulado no capítulo anterior, pois o vetor 
normal ao plano xOy é igual ao vetor k = (0,0,1) e uma curva neste plano tem vetor 
tangente igual a y'(t) = (a’(t), y'(6), 0). 

Portanto k x y'(t) = (—y' (t), x'(t), 0) e 


0 
VON) VANO) =| y) 2) 0 =0+y(07>0 
1 


Teorema 5.3.1 (Stokes) Sejam r(u,v) = (x(u, v), y(u, v), z(u,v)) a parametrização de 
uma superfície S regular e D =r(R) uma regiáo simples limitada por uma curva fechada 
regular por partes C orientada positivamente. 


Entao 
$ ras +Gay+Haz= | | roxas = ff rot(X).NdA 
C D D 


Apresentaremos a seguir uma demonstracáo do Teorema de Stokes. Antes porém, vejamos 
alguns exemplos e exercícios de aplicação. 


Exemplo 5.12 Se X(x,y,2) = (F(x,y),G(x,y),0) então rot(X) = (0,0,Gy(x,y)—Fy(x,y)) 
é ortogonal ao plano Oy. Seja C uma curva parametrizada por y(t) = (a(t), y(t),0) no 
plano xOy, orientada positivamente, que limita uma região R C R?. Pelo Teorema de 
Stokes, 


$ F(z,y)de + G(a, y)dy = / 3 rot(X)dS = / I [Gy (a(t), y(t)) — Fy(a(t), y(t) |. 


Este é precisamente o enunciado do Teorema de Green. Vemos assim que o Teorema de 
Stokes é uma generalização do Teorema de Green. 


Exemplo 5.13 fo ydx—axdy+ [cos(z*) — vyldz para C a intersecção do plano x+y+z = 1 
com o cilindro z? + y? = 1 orientada positivamente. 

Solução Usaremos o Teorema de Stokes. C limita a região R do plano definida por 
a? ty? <lez=1—a~—y (um gráfico). Um vetor normal ao plano é dado por (1,1,1) 
e o campo vetorial é igual a X = (y, —x,cos(22) — xy). Logo, rot(X) = (—2,y,—2) 
(verifique!). Portanto, pelo Teorema de Stokes, 


$ uds + —ady + [cos(2?) — xy]dz = [fics —2).(1,1,1)dA 


fvi + —ady + [cos(z”) — zydz = [fic + y — 2)dA. 
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Usando coordenadas polares, a região R é descrita pelas desigualdades 0 < 0 < 27 e 
0<r< 1. Temos assim, 


27 
$ ydx + —zdy + [cos(2?) — xy|dz = i / (—rcos(0) + rsen(0) — 2)rdrd0. 
C o Jo 


r3 3 1 


27 
f ydx + —zdy + [cos(2?) — zyļdz = / [—-—cos(0) + E sen(0) — 2r]| de = 
C 0 3 3 0 


27 1 1 
f ydx + —ady + [cos(2°) — zyļdz = I [-=cos(0) + = 
c 0 3 3 


Exemplo 5.14 Encontre fo x?ydz + y?zdy — 22adz para C a intersecção do plano x + 
y +z = 1 com os planos rOy, rOze yOz orientada pe 

Solução: Escrevendo o campo vetorial X(x,y,z) = (2?%y,y?*z,22x), temos V x X = 
(-y?, 2, —z"). Verifique que a curva C limita a região triangular D do plano z = 1-x—y 
definida pelas seguintes desigualdades: x > 0, y >O0e2z>0. C é um triângulo de vértices 
A = (1,0,0), B = (0,1,0) e C = (0,0,1). A orientação positiva em relação à superfície 
significa que o triângulo é percorrido de A para B, de B para C e de C para A. Verifique 
também que D é igual a imagem de triângulo R definido pelas desigualdades x > 0, y > 0 
ex+y<l. 

Utilizando o Teorema de Stokes e também que que z = 1 — x — y, temos: 


pe zdy — z an (V x X).NdS = 


= | [cv y)2, -LA 


Para obter os limites de integração, descrevemos a região R pelas seguintes desigualdades: 
0<x<1leQ0< 1— ze escrevemos: 


$ xr ydr + y?zdy — 22adz == I [-y? + (1-x- y? — a? ldA = 
C R 


= f [1-22 +2oyjda 
R 


Finalmente, calculando a integral dupla (a cargo do leitor): 


Exemplo 5.15 Seja X(x, y, 2) = (2%, —y?, x?) um campo vetorial. Encontre f fp rot(X)dS 
para D o hemisfério x? + y? + 2? =1, z > 0 com vetor normal apontando para cima. 
Solução: A região D é limitada pelo círculo C := La? + y? = 1,2 = 0}, que orientamos 
positivamente em relacáo ao interior da regiáo D. 
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Pelo Teorema de Stokes f [prot(X)dS = [o 2ºdx — y?dy + x7dz. 
Para calcular a integral curvilínea, usamos a parametrização 
6 ++ (cos(0), sen(0),0) do círculo para obter: 


27 3 6 
/ ¿dx — y dy + 1*dz = / [—sen*(0)cos(0)]d0 = sem ) an =0. 
C 0 
É interessante (mas mais longo) calcular diretamente f f prot(X)dS e confrontar os re- 
sultados. 
Demonstracáo: 


A apresentacáo desta demonstracáo tem como objetivo proporcionar um melhor enten- 
dimento dos conceitos e ideias envolvidos. Em particular, pretende-se enfatizar a inter- 
pretacáo de que estamos generalizando o Teorema Fundamental do Cálculo. Veremos 
como o Teorema de Stokes aplicado a uma situacáo especial, segue do Teorema de Green. 
Este por sua vez, é consequéncia do Teorema Fundamental do Cálculo, também aplicado 
em uma situacáo especial. A leitura desta demonstracáo náo é essencial, mas para utilizar 
corretamente o Teorema é preciso saber verificar as circunstáncias em que ele pode ser 
aplicado. 


Primeiramente desenvolvemos a integral de superfície do lado direito para esclarecer o que 
deve ser demonstrado: 


yl rot(X).NdA = TJ (Hy — Gz, F, — Hs, Gs — Fy) (tu X ro)dA 
D R 


Lembrando que 


rulu, v) = (£u lu, V), Yulu, V), zulu, v)) € rolu, v) = (zy (u, v), Yo (u, v), 210,0). 
então 


ru X Ty = (Yuzv — Yoru; Zuu — Zulu, TuYv — tn) 


O roteiro da demonstração é bastante simples. 
Calculamos cada uma das parcelas de 


Pact | Gay+ | Há: 
c c C 


Usamos o Teorema de Green em cada parcela para passar da integral curvilínea para a 
integral dupla. 

Finalmente, somamos as trés parcelas para obter exatamente o lado direito da igualdade. 
Faremos o cálculo para a primeira parcela, bo Fdx; para as demais, o cálculo é semelhante. 
Queremos mostrar que 


$ rao= | | trigo — zuv) + Felen — 22010 
O R 


Escreva a parametrizacáo de C da forma 


a(t) = r(u(t), v(t)) = (z(u), v(t)), yu), v(t)), z(u), v(t))). 


A integral curvilinea se escreve: 


$ Par = | FEO oO) vlo), (u(t), (O) (at 
C Q 
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Pela Regra da Cadeia 2'(t) = zuu’ (t) + xuv (t), portanto, substituindo temos: 


f Pac = f FEUD VO) E) vE), (ult), Ele (E) + 2.0" Olde 
Ou 


$ Rits fri E FeO aw Oa = f KaD t E a 


Esta última é uma integral curvilínea no plano (u,v) em um caminho a que limita uma 
região simples R. 
Nesta região, aplicamos o Teorema de Green: 


$ ral) Pedo | [Lt e) — EE e 0) molas 


Note que a primeira parcela é igual a 


o 
gu E (elu, v), y(u, v), z(u, v))£y| = [Fr Tu + FyYu + Fotu Lo + Fou 


Enquanto a segunda é igual a 


5 UF ((a(u, 0), y(u, 0), z(u, V))£u] = [E EE UF 2 lt + Ftw 


Subtraindo a segunda expressão da primeira, e omitindo o argumento (x(u, v), y(u, v), z(u, v)) 
de F, obtemos, após simplificações algébricas: 


o 


o 
ðu [F' xo] > = [Fra] = Fy Pulo = Duda] + E, [2 To = Tuzu] 


Ov 


Isto é: 
$ Fde = / fr, (tuto — uty) + El tuto — Luto) dA 


Procedendo analogamente para as demais parcelas temos: 


$ Gdy = 1 [ [Gx (uy — ToYu) + Gz (ZuYo — 2vYu)|dA 


$ Hdz = Ir el Duty = Luu) + Hy(Yuzv Yu) OA 


Finalmente, somando as três parcelas e coletando os termos comuns obtemos 


f Fdz + Gdy + Hdz = 
C 


= | | Fs- Gina Ue [Fe Hal [Sha Zola] | [en Fy | (tut — Lo A 


Concluimos assim que 


$ Pac + Gdy+ Haz | | rot(x).waa, 
C D 


Como queríamos demonstrar. 
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5.4 Campos conservativos no espaço 


Vimos na Seção 4.7, Teorema 4.33, que a Integral Curvílinea de alguns campos no plano ou 
no espaço não dependem do caminho y(t), a < t < b ao longo do qual a integral é calculada 
mas apenas dos pontos inicial e final y(a) e y(b). Estes campos, chamados conservativos, 
são importantes na Física pois o trabalho independe da trajetória. Vamo a seguir. utilizar 
o Teorema de Stokes para caracterizar os Campos Conservativos no Espaço. 


Teorema 5.16 Seja X = (F,G, H) um campo vetorial no espaço com componentes contínuas 
e diferenciáveis em uma aberto U C RÌ, exceto possivelmente em um número finito de pon- 
tos. Então são equivalentes: 


1. bo X.dr = 0 para toda curva fechada C regular por partes contida em U. 


2. A integral curvilínea de X do ponto À até o ponto B independe da curva regular por 
partes que liga A a B. 


3. X é um campo gradiente de alguma função f em U. 


4. rot(X) = 0 


Demonstração: A demonstração deste teorema é praticamente idêntica a demonstração 
do Teorema 4.7.1. A maior diferença é que aqui utilizamos o Teorema de Stokes em vez 
do Teorema de Green. Vamos mostrar que (1) => (2) => 3 => (4) => (1). Inicialmente 
para ver que (1) => (2) , considere duas curvas C4 e C2 regulares por partes contidas em 
U e unindo os pontos A e B (ver figura). 


Pela Observação 4.23, dada a curva Ca, existe a curva —C uma curva idêntica a C2 mas 
com ponto inicial B e ponto final A. Então C = C1 U —C> é uma curva fechada e por (1) 
temos que: 


$ Xr = X.dr — X.dr = 0. 
C Cı Ca 


Segue que e, X.dr = fe, X.dr. 

Para mostrar que (2) => (3), utilizamos um argumento análogo á demonstracáo do item 
correspondente no Teorema 4.7.1, que omitiremos aqui. 

A implicacáo (3) => (4), será vista no Exercício 5. 

Finalmente vamos mostrar que (4) => (1). Para isto seja C uma curva fechada regular 
por partes e seja S uma superfície cuja fronteira é C. Vamos supor que se X possui pontos 
excepcionais onde o campo náo é contínuo ou diferenciável que S evitas estes pontos. 
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E sempre possivel conseguir isto, mas nao demonstraremos este fato aqui. Pelo Teorema 


de Stokes temos 
$ Xar= ff rot(x) a5 = 0. 
C S 


Exemplo 5.17 Calcule a integral do campo X (x,y,z) = (yz +x, £z +y, ry +2?) ao longo 
da curva C obtida como interseção da semiesfera £? + y? + 2? = 1,2 > 0 com o plano 
y=0. 

Solução: A curva C é um semicírculo no plano y = 0. Consideremos a curva orientada 
de (—1,0,0) para (1,0,0). Observe que rot(X) = 0 pois 


i j k 
ro(X)=V xX = E 5 E =(x—2,—-(y—y),z— 2) =0. 
Y2+txu az+Hy xy + 2” 


Pelo teorema anterior existe uma função potencial f para o campo X, ou seja, tal que 
Vi = es a, of) = X. Vamos calculá-la. 


Como of = yz + x temos que f(x,y,z) = yza + = + A(y,z) onde A(y, z) é uma função 


que não depende de x. Analogamente trabalhando com as outras derivadas parciais de f, 
2 

que conhecemos, obtemos: f(x,y,z) = yzz +% +A(x, z) e f(x,y,2) = yzz + = + A(x,y). 

Comparando vemos que uma função potencial para X pode ser dada por: 


2 2 3 


f(x, u,2) = yzz +4 5 5 + a 


Como também pelo teorema anterior sabemos que a integral não depende do caminho, 
mas apenas dos pontos inicial e final temos que: 


f X.dr = f(—1,0,0) — f(1,0,0) =0. 
C 
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5.5 


EXERCICIOS 


. Encontre o fluxo do campo X (z, y, z) = k sobre a esfera x? + y? + 22 = 1, com vetor 


normal apontando para o exterior da bola. 


. Encontre o fluxo do campo X (x,y,z) = —xi—yj—zk sobre o cubo [0, 1] x [0,1] x [0,1] 


com orientação para fora. 


. Encontre o fluxo do campo X(x,y,2) = (xy,yz,2x) sobre a parte do gráfico de 


z = 4 — 2? — y? que está sobre o quadrado 0 < x < 1e0< y< 1, com orientação 
para cima. 


. Para cada campo de vetores abaixo, encontre rot(X ): 


(a) X(7,9,2) = (x, y, 2) 
(b) X(x,y,2) = (xy, yz, zx) 
(O Xr yz) ere ne) 


. Seja X(x,y,2) = Vf(x,y,z) o campo gradiente para f(x,y,z) uma função com 


derivadas parciais de ordem dois contínuas. Mostre que rot(X) = 0. 


. Utilize o teorema de Stokes para calcular: 


(a) 
i ydx + zdy + zdz 
CG 
onde C é a curva de interseção do plano z +y = 4 com a a esfera x? +y? + 2? = 
A(x +y). 
(b) 
f ydx + zdy + zdz 
C 
onde C é a curva de interseção do plano z +y = 2 com a a esfera x? +y? + 2? = 


2(x + y). (Sugestão: faça um esboço da curva de interseção e da área que ela 
envolve e observe que esta área é facil de calcular.) 


. Encontre f g€ *dx +e*dy + e*dz para C a interseção do plano 2x + y + 2z = 2 com 


os planos coordenados, orientada positivamente. 


. Encontre [.(y— z)dz — (z—zx)dy+ (x-—y)dz para C a interseç ao do plano +2 = 1 


com o cilindro z? + y? = 1 orientada positivamente. 


. Seja F(x,y, z) = (2zz + y?, dry + Sy, e” + 27). 


(a) Verifique se F é um campo conservativo. 
(b) Determine uma função potencial para F. 


(c) Seja C a curva obtida como interseção da superficie z = 9 — x? — y?, z > —4 
com o plano y = 2. Calcule le F.dr especificando uma orientação para C. 
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DEST 


Teorema da Divergéncia 
(Gauss) 


AULA 6: TEOREMA DA DIVERGENCIA (GAUSS) 


Objetivos 6.1 Os objetivos desta Aula sao: 


e introduzir os conceito de divergéncia de uma campo vetorial: a divergéncia é uma 
taxa de espalhamento de um campo vetorial; 


e relacionar o fluxo de um campo vetorial através de uma superficie fechada com a 
integral tripla da divergência do campo na região limitada pela superfície (Teorema 
da Divergência ou Teorema de Gauss); 


6.1 Divergéncia de um campo 


Nesta aula será exposta a segunda generalização do teorema de Green que foi proposta no 
início do capítulo anterior. Isto é, a que relaciona o fluxo de um campo sobre uma superfície 
S com uma integral tripla no interior da região do espaço limitada pela superfície. 


Voltemos à situação bi-dimensional: imagine que temos um campo vetorial X = (F,G 
em R? e uma curva simples e fechada, C, imagem de uma parametrização regular y(t) = 
(x(t), y(t)). 


O Teorema de Green nos diz que 


f| -Gas+Fay= | f |F +GoldA. 
C R 


Observe o sinal e a ordem em que foi escrita a integral curvilínea à esquerda da igualdade. 
Esta integral curvilínea pode ser interpretada do seguinte modo: Se y(t) = (a(t), y(t)), 
então y (t) = (a’(t), y/(t)) é o vetor tangente ao caminho e o vetor (y'(t), —x'(t)) é vetor 
ortogonal ao caminho. Além disso, o produto escalar do campo X com o vetor normal 
é igual a X.n = Fy’ — Gz’. O sinal do produto escalar X.n nos diz se o campo está 
“entrando” (negativo) ou “saindo” da região R. Portanto, a integral curvilínea Lo —Gda + 
Fdy é o fluxo do campo sobre a curva. O fluxo pode ser interpretado como uma medida 
do espalhamento do campo em relacáo a curva. Isto é de certa forma medido pela funcáo 
F; + Gy que aparece na integral dupla à direita. 

As considerações acima motivam a seguinte definição: 


Definição 6.1.1 Chama-se divergência de um campo vetorial X = (F, G, H) a seguinte 
função escalar: div(X) = Fr + Gy + H; 


Exemplo 6.2 Se X(x, y, z) = (x,y,z), então div(X) = F, + Gy + H; =1+1+1=3 


Exercício 6.3 Se X (x,y,z) = (—x, —y, —z), calcule div(X) 
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6.2 Teorema da Divergência de Gauss 


Estamos prontos para enunciar Teorema da Divergéncia , também conhecido como Teo- 
rema de Gauss, que relaciona o fluxo de um campo sobre uma superficie com a divergéncia 
na regiao limitada pela superficie. 


Teorema 6.2.1 (Divergéncia de Gauss) Seja X (x,y,z) um campo vetorial com deri- 
vadas parciais contínuas definido em uma região U C R*. Suponha que Q seja uma região 
contida em U limitada pela união finita de superfícies regulares orientadas com vetor nor- 
mal apontando para fora da região Q. Então, 


I [ X.NdA = / / A div(X)aV 


Antes de ver uma demonstração do Teorema da Divergência, vejamos alguns exemplos e 
exercícios. 


Exemplo 6.4 Considere a superficie S = S1 U S2 U S3 regular por partes, definida por: 
Sı = {£2 +y? =1;0 < 2 < 1}, S2 := {£2 +y? < lie = 0} e S3 :=l0*+y4<1i2=1). 


A orientação de S ” para fora” é obtida tomando-se os seguintes vetores normais: 

N = (—sen(0), cos(0),0) em S1; No = (0,0, —1) em S2 e N3 = (0,0,1) em S3. 

Seja X(x,y,2) = (x,y,z) ( campo radial), com div(X) = 1+1+1 = 3. Pelo Teorema 
da Divergência o cálculo do fluxo de X sobre S é igual à integral tripla da divergência no 
cilindro sólido 2 limitado por S, isto é: 


f [xxi f | [atoa = f f f av = vom) = 


Exemplo 6.5 Sejam y(x,y,z) = (x7, y?, z?) e S = Sı U S2 U S3 a mesma superfície do 
exemplo anterior. 
Como div(X) = 2x + 2y + 2z, pelo Teorema da Divergéncia: 


JIGXNdA= f f Jodiv(X \dV =2[] fala + y + 2Jav 
Usando coordenadas cilíndricas, x = rcos(0), y = rsen(0), z 


27 
uid x+y + zļdV = 2/ i. fi rcos(0) + rsen(0) + z|rdzdrdo 


= 1 
= =2 f fr (cos(0) + sen(0)) + slárdo = = re [3 (cos(@) + sen(0)) + alarde = =” 
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Exemplo 6.6 Encontre o fluxo do campo X(x, y, 2) = (cos(y) + z°)i+ (sen(az) + y?)j + 
(xe*Y + z)k sobre a superfície S, regular por partes, definida pelas equações x = 4 — y?, 
xz + z = 6, o plano yOz e rOy com vetor normal apontando para fora da superfície. 
Solução: Faça um esboço da superfície S. Lembre-se que uma superfície em R? definida 
por uma equação que envolve apenas duas variáveis é uma superfície cilíndrica. Deste 
modo, por exemplo, a equação x = 4 — y? representa um cilindro parabólico. 

Estamos em condições de usar o Teorema da Divergência, pois o campo X está definido 
em R3. Como div(X) = 2y+1, se denotamos por £ a região interior à superfície S, então: 


f fxxas= | | [Cu+nav 


Q pode ser descrito pelas desigualdades: 0<z<6-=x,0<xw<4-y?e-2<y<2 
(verifique isto!) 
portanto, 


2 py? p6-x 
T X.NdS = / / (2y + 1)dzdzdy. 
S -2 Jo 0 
4—y? 6-2 4—y? 
L / / (2y + 1)dzdxdy = F / [(2y + 1)(6 — x)|drdy 


4— 
L. f É [(2y + 1)(6 - a)jdzay = f" (2y + 1)16(4— ya) — E É Jay 


: 2 (4 — yy 2 5 y? 3 2 
(2y + 1)[6(4 — y”) > ldy= | [-y? - 5 — Ay? — 24? + 32y + 16]dy 
-2 3 
Observe que o intervalo de integração é simétrico, |—2, 2], portanto a integral das funções 
ímpares é nula. 
Logo 


2 4 
694 
f [-y® — © — ay — 2y? + 32y + 16]dy = —— 
= 2 15 


Exercício 6.7 Encontre o fluxo do campo X (x,y,z) = (x,y,z) sobre a esfera com centro 
(0,0,0) e raio igual a R, com normal apontando para fora. 


Demonstração do Teorema da Divergência: 

Analogamente ao que foi feito na demonstração do Teorema de Green, vamos supor que 
que a região Q no espaço seja simples e que possa ser escrita simultaneamente como uma 
região entre gráficos, entre em cada uma das direções Ox, Oy ou Oz. 

Isto significa que Q se expressa, simultaneamente da seguinte forma: 

Q :=([h1(y, 2) < x < gi[y, z), (y, 2) € Ri CR?) com R; a região do plano (y, z) ou z = 0 
sobre a qual projeta-se a região sólida Q ou 

Q:=(ho(x,2)<y< g2(z, z), (2,2) € Ra C R?], Ro imagem de Q pela projetção ortogonal 
sobre o plano de coordenadas (x, z) ou 

Q := [h3(x,y) < z < ga(x, y), (x,y) € R3 C R?}, R3 imagem de Q pela projetção ortogonal 
sobre o plano de coordenadas (x,y). 


AULA 6: TEOREMA DA DIVERGÉNCIA (GAUSS) 


93 


94 


A demonstração neste caso é similar à que foi apresentada para o Teorema de Stokes. 
Isto é, se X = (F,G, H) = Fi+Gj+ Fk, calculamos cada uma das parcelas de f f fo[Fr+ 
Gy + H.|dV separadamente e provamos que 


f| [r= | [ rinda 
fJ [ewv=f [eines 
J| [nae | [tunas 


Em cada um dos casos usamos a descrição mais conveniente de Q, para que possamos usar 
o Teorema Fundamental do Cálculo. 

Vemos também como é importante o Teorema de Fubini, que permite-nos escolher a ordem 
de integração. 

O resultado é obtido somando-se as três igualdades acima. 

Provemos a primeira igualdade, as outras seguem de maneira análoga: 

Escrevemos dA, para destacar que a ordem em que calculamos a integral repetida: pri- 
meiramente calculamos a integral em relação a x, em seguida calculamos a integral dupla 
sobre a região Rı contida no plano (y, z). 

Ela indica também que a integral de superfície é calculada sobre a imagem das superfícies 
parametrizadas (y, 2) > (ha(y,2),y,2) e (y, 2) > (g1(y, 2), y, 2) respectivamente. 


gı (y,z 
[[fev-ff ] aid 
Ri Shi(y,z) 


Pelo Teorema Fundamental do Cálculo: 


J | [eave f f Eo: -Etv 24. 


Por outro lado, o vetor unitário da normal N aos gráficos de gı (y, z) e de hı (y, z) são os 
ty et O, ð Ohi Oh 

vetores unitários de ( (ppt de) e de (—1, SH, 2). 

Por que a diferença dos sinais? Porque estamos tomando o vetor normal apontando para 

fora da região 2. 


Assim a ea para o produto escalar Fi.N sobre cada um dos gráficos é: 
F(gi(y, 2), y, 2) e ~F (hı (y, 2), Y, 2) 


/ [ FiNdA = f [E [Pau 2) y2) = Fla), 2)]d Aya 


Usando a igualdade, temos 


Pf Ew f f Eon -Eyd f f Finas 


Como foi dito anteriormente, o mesmo raciocínio é empregado para obter as demais igual- 
dades. Concluindo a prova do Teorema da Divergência para o caso especial da região Q 
acima. 
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Para o caso de regiões mais gerais, seguimos o que foi feito no teorema de Green, por meio 
da decomposição da região dada em regiões do tipo acima. 
Para isso, a orientação exerce um papel fundamental. 


6.3 Teorema da Divergência de Gauss em regiões mais ge- 
rais 


Considere o seguinte exemplo. 


Exemplo 6.8 Vamos tentar aplicar o Teorema 6.2.1 para calcular o fluxo do campo 
X(x,y,2) = (= a ise ri e DE o =z) sobre a esfera centrada na origem e de raio 
2 com vetor normal apontando para fora. Ao fazer isto temos um problema já que o 
campo X não está definido na origem e portanto o Teorema 6.2.1 não se aplica a este 


campo, se nos considerarmos que (2 é a esfera de raio 2. 


Podemos utilizar o Teorema da Divergência de Gauss, em regiões mais gerais do que as 
que consideramos até agora. 
Seja Q a região do R? definida por 


Q ={(z,y,z) E€ R tal que 1 <a? +y? +2? <4 


Observe que a fronteira 90 de 2 é agora constituida por duas esferas uma exterior S de 
raio 2 e uma interior Si de raio 1. Como o campo X só não está definido na origem e 
satisfaz as hipóteses do Teorema 6.2.1 podemos aplicar o Teorema de Gauss nesta região 
mais geral. 

Para isto vamos utilizar o Exercício 5.10 que mostra que o fluxo do campo X é 47 para 
qualquer esfera de raio a independentemente do raio da esfera. Faça agora o seguinte 
exercício: 


Exercício 6.9 Mostre que a divergência do campo X na região Q acima é zero. 


f fax nia= [ [o xxaa= [| f diwa =o 


ai x.Nda= f x.vaa+ f X.NdA 
S1US2 Sy Sa 


ef fs, X.NdA = — [ Ss, X.NdA = 47, um resultado que já conhecíamos. 


Temos então: 


Mas como 
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Mas agora podemos obter um resultado bastante surpreendente. Seja (2 a região do espaco 
limitada internamente pela esfera Sı de raio 1 e externamente por qualquer superfície S 
tal que Sı está no interior da região limitada por S. 

Podemos facilmente utilizando o que já fizemos calcular f f g X.NdA. Pelo Teorema da 
Divergência de Gauss temos: 


/ X.NdA= / X.NdA= / / | dnl av =0 
90 SUS Q 
I x.Naa= f xNda+ | [ x.vaa 
S1US S1 S 


e ÍS Ís, X.NdA = —47, Temos 
J | XNis=4r 
S 


Mas como 
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6.4 EXERCICIOS 


1. (a) Se X(x,y,2) = (2, y, x), mostre que div(X) = 0 
(b) Se X(x, y, z) = (—y, 2,0) mostre que div(X) = 0. 


2. Vamos utilizar a notação V = (2, = 2) 


(a) Prove que se U é uma funcáo com derivadas de segunda ordem contínuas, entáo 
rot(VU) = 0. 

(b) É verdade que div(VU) = 0? 

(c) Mostre que podemos escrever div(X) = V.X 

(d) Se X(x,y, z) = ( 


2x 
r? +y? +z 


7) 72 Td ae a =), calcule div(X). 
3. Prove que div(rotX) = 0 


4. Calcule o fluxo do campo X (x,y, z) = 1y?*j+y2?+2x*k sobre a esfera x?+y? +2? = 1, 
com vetor normal apontando para o exterior da bola. 


5. Use o Teorema da Divergéncia para calcular o fluxo do campo X (z, y, z) = (1*,0, zx?) 
sobre a esfera centrada na origem e de raio 1 com vetor normal apontando para fora. 


ee 


E 7, = 2x 2y 2z 

(a) E possível calcular o fluxo do campo X (z, y, 2) = (= a y a 3) 
sobre a esfera centrada na origem e de raio 1 com vetor normal apontando para 
fora, utilizando o Teorema de Gauss? Por que? 


(b) E na esfera (x — 2)? + y? + 2? = 1? Explique. 


(c) É verdade que os fluxos do campo X nas esferas a24+y?+2? = lea?+y?+22 = 4 
sáo iguais? Explique. 
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